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Course information 



Course outline 

Aim: to introduce 1st year physics graduate students to the theory of 
quantum fields. 

• Relativistic wave equations (Ryder Ch 2, G & R Ch 4.1, 6.2). Rela- 
tivistic notation; Klein-Gordon equation. 

• Lagrangian formulation (Ryder Ch 3, G & R Ch 2). Lagragian particle 
mechanics; Euler-Lagrange for real scalar field; Noether's theorem for 
scalar fields. 

• Canonical quantisation (Ryder Ch 4, G & R Ch 4/7). Canonical com- 
mutation relations for real scalar field; Fock space; Number, energy and 
momentum operators; Propagators. 

• Fermions (Ryder 4.3, 6.7). Dirac equation; Lagrangian formulation; 
Canonical quantisation; 

• Interacting fields 1 Interaction picture; S-Matrix and transition ampli- 
tudes; Wick's theorem; Matrix elements in perturbation theory; Scat- 
tering cross section 2^2 particles. 

• Scalar Quantum Electrodynamics Maxwell equations; Lagrangian for- 
mulation of Electromagnetic field equations; Canonical quantisation of 
EM field in Lorentz gauge; Gauge symmetry; Feynman rules for SQED; 
photon polarisation sums in scattering probabilities. 

• Scalar field theory in an expanding Universe Friedmann-Robertson- 
Walker metric; Lagrangian; 

Learning outcomes 

By the end of the course the student should be able to: perform four- 
vector algebra; derive the Euler-Lagrange equations for any field theory; 
understand and apply Noether's Theorem; give an account of the canonical 
quantisation procedure for fields; explain how transition amplitudes are cal- 
culated in field theory; write down and calculate simple Feynman graphs; 
give an account of quantum fiuctuations of the field in the vacuum. 
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Teaching methods 

There will be 3 lectures a week throughout the term. Problem sheets will 
be given every two weeks. 
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Assessment 

MSc students Problem sheets will count for 40% of the total mark for the 
course, the remaining 60% will come from the 1.5 hour exam at the 
end of the year, which will have a choice of two questions from three. 

DPhil students Problem sheets will count for 100% of the total mark. 
Reading list 



Field Quantization, W. Greiner and J. Reinhardt (Springer, Berlin, 
1996), 

Quantum Field Theory, F. Mandl and G. Shaw (Wiley, Chichester, 
1984). 

Quantum Field Theory, L. Ryder (C.U.P., Cambridge, 1984). 

An introduction to quantum field theory, Michael E. Peskin, Daniel V. 
Schroeder (Addison- Wesley, Reading, Mass; 1995). 

Introduction to Gauge Field Theory, D. Bailin & A. Love (Adam Hilger, 
1993) 

Quantum field theory: a modern introduction, Michio Kaku (O.U.P., 
Oxford, 1993) 

Fields, W. Siegel (http://insti.physics.sunysb.edu/ siegel/plan.html). 

The books are complementary in some way. Greiner & Reinhardt is 
extremely thorough. Ryder has a good balance of topics and starts fairly 
simply, but used path integral quantisation which is not the most pedagogical 
way to introduce field theory. Mandl and Shaw is pedagogical but (like 
Greiner & Reinhardt) adopts different Fourier transform and normalisation 
conventions which can confuse the unwary. Peskin & Schroeder covers the 
most material and is probably the book of choice for serious students of field 
theory. There are quite a few errors: a list can be found at 
http : / /www . slac . Stanford . edu/~mpeskin/QFT . html 
Bailin & Love starts at a higher level than this course, but is concise 
and has a unique discussion of the S-matrix generating functional. Kaku is 
very complete but rather rushed. Siegel's PDF lecture notes are the most 
complete exposition I have seen. They are very advanced and will cost you 
a lot to print at over 800 pages. 
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Prerequisites 

The course assumes that you have aheady encountered 

• Quantum Mechanics (Schrodinger equations, free particle solution, har- 
monic oscillator). 

• Special Relativity (Lorentz transformations, space-time interval, Lorentz 
invariance) . 

• Electromagnetism (Maxwell's equations in free space). 

• Cauchy's theorem (contour integration). 

• Fourier transforms. 

• The Dirac 5-function. 

• Scattering cross-sections. 

It will help to have come across the following subjects, mostly covered in the 
Sussex undergraduate course Theoretical Physics II, as only a brief review is 
given here. 

• Lagrangian formulation of classical mechanics (e.g. Classical mechan- 
ics, T.W.B. Kibble and F.H. Berkshire, London: Longman, 1996). 

• 4-vector notation in Special Relativity (e.g. Chapter 2 in General rel- 
ativity and cosmology, T.L. Chow, Winnipeg: Wuerz, 1994). 

• Natural Units (Mandl and Shaw, Section 6.1). 



Course Lecturer 

Mark Hindmarsh, Scitech, Arundel 213. Telephone: 8934 
E-mail: m . b . hindmarsh [NOSPAM] Sussex .ac.uk 
Course web: http : //www . pact . cpes . Sussex . ac . uk/~markh/RQFl/ 
Office hour: Tuesday 2-3pm, or by arrangement (email is best). 
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1 Preliminaries 



The reason for studying quantum field theory is that it is the theory which 
best describes elementary particles and their interactions. Field theory rec- 
onciles quantum mechanics and special relativity and is automatically a 
quantum-mechanical theory of many particles. Where we can calculate with 
great precision, such as in Quantum Electrodynamics, spectacularly accurate 
agreement with experiment is found. However, field theory is not entirely 
satisfactory. Compared to the elegant geometrical framework of General 
Relativity, it is a mess, and when one first tries to use it to calculate simple 
physical quantities, we immediately encounter infinite results. Making sense 
of these infinities occupies a large proportion of any quantum field theory 
course, but fortunately most of them can be ignored if we neglect interac- 
tions between particles. There are still non-trivial physical predictions to 
be made in non-interacting (free) field theory, associated with non-trivial 
boundary conditions, such as in the Casimir effect, or with strong gravita- 
tional backgrounds, such as we expect to find near black holes or in the early 
Universe. 

This course is aimed towards computing the amplitude of fluctuations in 
a real scalar field in an exponentially expanding space-time, and does just 
enough (but no more) to get there. We will prove the famous result that in 
this spacetime, which is a limiting case in the class of inflationary universes, 
the power spectrum of the fluctuations in a massless scalar field is (iJ/27r)^, 
where H is the expansion rate. The scalar field which causes inflation, the 
infiaton, is effectively massless during inflation, and these vacuum fluctua- 
tions form the basis of the explanation of the perturbations we observe in 
the Cosmic Microwave Background today. 

1.1 Special Relativity 

We start by quoting a formulation of the Principle of Special Relativity: 

Fundamental physical laws are the same for all observers moving 
with constant velocity relative to one another. 

Special Relativity is a rather peculiar theory as it appears to be a theory 
about other theories. However, when applied to dynamics it does have 
testable consequences, which have of course corroborated the theory to great 
accuracy. It is in fact true in only a restricted domain, where gravitational 
fields are weak and can be neglected. Although it may not seem like it, 
the gravitational field at the surface of the earth it sufficiently weak to be 
neglected for most purposes in relativity. 
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1.2 Lorentz transformations and space-time interval 



Suppose a frame of reference F' is moving at velocity v in the x direction 
relative to another F. Special Relativity tells us that coordinates measured 
in F' are related to coordinates measured in F by a Lorentz transformation 

t^t' = -f{t-vx/c^), 
x^x' ^ ^{x-vt), 

y ^y' = y, 

z ^ z' — z, 

where 

7 = 1/^1 - v^lc". (1.2) 
Although neither the distance nor the time interval between two events is 
observer-independent in special relativity, there is a concept of space-time 
interval, which contains a little of both, and is something that all observers 
can agree upon. In its infinitesimal form, the interval ds between two events 
at (t, x, 7/, z) and {t + dt,x + dx, y + dy, z + dz) is given by 

ds'^ = c^dt^ - dx^ - dy"^ - dz^. (1.3) 

One can easily check that the value of ds docs not change under this trans- 
formation. We say that ds is Lorentz invariant. 



1.3 Relativist ic notation 

If a theory - which is usually formulated in terms of a set of equations 
- is consistent with the Principle of Special Relativity we often say that 
it is Lorentz covariant. There is a notation which, if followed correctly, 
automatically ensures that equations are Lorentz covariant. Quantities are 
assembled into ^-vectors (and later on we will encounter ^-tensors) which 
transform in a simple linear way when one compares their values between 
observers in different states of uniform motion. Just as one can form a 3- 
component vector from three spatial coordinates 

x" = (a;\ a;^, x^) = {x, y, z), (1.4) 

one can define the 4-vector space-time coordinate for an event by 

x^ = {x^ , x^ , x'^ , x"^) = [ct,x,y,z). (1-5) 

We shall adopt the following conventions when labelling vectors: 

• Greek letters from mid-alphabet for space-time indices (/x, z/, p, cr, . . .); 

• Roman letters from mid-alphabet for spatial indices (i, j. A;, . . .); 

• bold face will also be used for 3-vectors, e.g. x. 
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1.4 Matrix representation of Lorentz transformations 



Recall the form of a Lorentz transformation between two frames of reference 
moving with velocity v in the x direction relative to one another (1.1). This 
linear transformation can be represented by matrices acting linearly on 
the coordinates: 



X' 



or in differential form 



u=0 



(1.6) 



i/=0 



The matrix A has entries 



A^- 



/ 7 



V 



-■jv/c 



-'fv/c 7 



1 • 



1 / 



;i.7) 



;i.8) 



where fj, labels the rows and u the columns. 

In order to save writing a large number of summation signs, we often use 
Einstein's convention that repeated indices in an expression are summed over 
automatically without the need for a summation sign. In this convention, 
Equation (1.7 becomes 



dx" dx'^ = A'*,dx^ 



(1.9) 



1.5 Space-time metric 



It is convenient in special relativity (and fundamental in general relativity!) 
to define a matrix gf^^,, which is used in the expression for the space-time 
interval: 



ds'^ = ^ gf.^.dx^dx". 

;t,f =0 



(1.10) 



We call Qf^i, the metric. In Special Relativity, the metric always takes a 
particular constant value. 



1 



(1.11) 



-1 / 
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We often call 77^1^ the Minkowski metric. In Einstein's repeated index con- 
vemtion, (1.10) can be written as 

ds'^ = rjf.^dx'^dx'' . (1.12) 

Note that the position of the indices is important. There is a difference 
between vectors with superscript indices and those with subscript indices, 
which shows up in their properties under Lorentz transformations. We may 
use the metric tensor to define an infinitesimal coordinate with a lower index: 

dx^ = rj^udx" (1-13) 

Vectors with subscript indices are termed covariant, to distinguish them from 
their contravariant partners with superscript indices. They transform oppo- 
sitely to contravariant 4- vectors, with the inverse of the Lorentz transforma- 
tion matrix A: 

dx^ dx'^ = {A-^)\dx^. (1.14) 

One can check that this is true by noting that the space-time interval can be 
written ds"^ = dx^dx^, and then substituting the transformation laws (1.7) 
and (1.13). One finds 

ds^ ds'^ = dx'^dx'^ - {K-^y ^dx^M pdxP . (1.15) 

Note that we have not just blindly substituted (1.7): we have the repeated 
index from a z/ to a p. The meaning is still the same: that index is to be 
summed over the values 0,1,2,3. However, if we have left the index as v 
equation (1.15) would be ambiguous, as we would not know how to pair off 
indices in the summations. This is an important rule with index notation: 
never use repeated indices twice on the same side of an equation. 

Continuing with equation (1.15), we note that the indices are paired 
and can be summed over. We can see that we are multiplying a matrix A 
by its inverse A^^ and therefore must obtain the identity, which expressed in 
index notation is 

{K-')\A% = (1.16) 
Here we introduce the Kronecker delta, defined by 




Note that the spatial components of a 4-vector with its index lowered (or 
a covariant vector) have the opposite sign to its counterpart with a raised 
index (a contravariant vector): 

Xq = a;°, Xi = -x\ (1-18) 
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Lastly, one can also define the metric tensor with raised indices as the matrix 
inverse of the covariant metric tensor 

^^" = 5^, (1-19) 

1.6 General Lorentz transformations: the Lorentz group 

One can explicitly verify that the transformation law (1.1) leaves the space- 
time interval ds invariant. By choosing space coordinates so that the relative 
velocity of two inertial frames is along the x direction, it follows that all 
Lorentz transformations leave the interval invariant. Let us see what we can 
infer about the matrices A from this condition. The interval transforms as 

ds"^ ds'^ = r^^^A^A'^c^x^dx^ = ds"^. (1.20) 

But ds^ may be written ds'^ — r]p„dx'^dx'^: thus we may infer that 

ll,uKK = Vpa. (1.21) 

Hence, any matrix A which leaves the metric 7/^,^ invariant under the trans- 
formation (1.21) represents a Lorentz transformation. These matrices form a 
group of transformations known as the Lorentz group. When combined with 
translation symmetry, x^ — > x'^ — x^^ + a^, with a constant 4-vector, it 
forms a larger group known as the Poincare group. 

1.7 Derivatives 

It is through studying derivatives that we discover the utility of the idea of 
covariant 4- vectors. Coordinate vectors are naturally contravariant, and it 
seems slightly perverse to introduce a covariant version with lowered indices. 
However, it turns out that there are 4- vectors which are naturally covariant, 
and the principal among these is the partial derivative 4-vector. 

The standard partial derivatives with respect to spatial coordinates can 
be collected together into a 3-vector: 

^ d ( d d d 
' dx'^ \dx^ ' dx"^ ' dx'^ 

The di are the components of the gradient operator V in an orthonomal 
basis e*, or 

V = e'di. (1.23) 



(1.22) 
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The relativistic generalisation is a covariant 4- vector 9^, which is defined as 

The transformation law for may be found from the chain rule: 

d d dx^ d 

'dxt^ ^ 'dx^ ~ dx^dx^ 

However, by differentiating 



;i.25) 



which follows from equation (1.6), we find that 



x'' = (A-^Y x"", (1.26) 



Hence the 4-vector differential really is a covariant 4-vector. 
1.8 4- vectors and the scalar product 

To recap, a contravariant 4-vector can be defined as any collection of four 
quantities which transform the same way as dx^ under a Lorentz transfor- 
mation. Similarly, a covariant 4-vector can be defined as a collection of 
four quantities which transform the same way as d/^. One can map any con- 
travariant vector into a covariant one with the metric 77^,^, an operation called 
lowering the indices, or one can go in the opposite direction with the inverse 
metric 1]'^'' (raising the indices). 

There are many important quantities than can be assembled into 4- 
vectors: for example, energy E and momentum p of a particle belong to- 
gether in a single momentum 4-vector p^, defined by 

p'^ = (£;/c,p). (1.28) 

We are used to the interval ds"^ = dx^dx^ being a Lorentz invariant quantity, 
but it is straightforward to see that one can form a Lorentz invariant quantity 
from any pair of 4- vectors a^, 6^. This is the scalar product, written 

a • 6 = / = ai'V'r]^^ = a^U" = - a^b^ - - (1.29) 

We can take the scalar product of the momentum 4-vector with itself and 
recover the relativistic relation between energy, momentum and mass: 

p2 = £;Vc2-pV = mV, (1.30) 
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Another invariant can be constructed by contracting the momentum 4- vector 
with the position 4-vector 



p - X = Et — p'X. 



(1.31) 



Another important object is the velocity 4-vector. In Newtonian mechanics, 
the velocity is dx^/dt, so it might be thought that dx'^/dt is the analogous 
quantity. However, dt is not Lorentz invariant, so the object as a whole does 
not transform like dif^ as a 4-vector should. Instead, we should differentiate 
with respect to a Lorentz invariant quantity. There is one ready-made for 
us in the form of the space-time interval ds, from which we can define the 
proper time dr — ds/c. The 4- velocity of a particle can then be written as 



1.9 The d'Alembertian 

There is a second order differential operator which students of electromag- 
nctism will already have come across, called the wave operator or the d'Alembertian. 
It is sometimes given its own symbol □: 



This is already Lorentz invariant, for it may also be written as d^j^d^rj^^ . You 
may recall that it was partly consideration of the theory of electromagnetic 
waves that led Einstein to formulate the special theory of relativity: elec- 
tromagnetism (in free space) is automatically a theory which is consistent 
with the special theory of relativity, a fact which becomes obvious when it is 
written down in terms of 4- vectors and tensors. This is not at all obvious if 
it is written down as Maxwell originally did, component by component. 

1.10 Lorentz covariance 

The power of 4-vector formalism is that when one uses it consistently to 
write down equations, the equation will keep the same form under a Lorentz 
transformation: that is, it will be Lorentz covariant. For example, we might 
want a relativistic generalisation of Newton's Second Law, F = dp/dt. We 
have already seen that differentiations with respect to time should be changed 
to differentiations with respect to proper time, so the analogous equation in 
relativistic dynamics should be 



□ 





dp^ 
d7 



(1.34) 



14 



The zeroth component of the force 4- vector, F°, is the (proper) rate of change 
of the energy of the particle. 

A potential source of confusion is whether physical quantities are rep- 
resented by covariant or contravariant vectors. As a rule, coordinates and 
momenta naturally have their indices up and derivatives naturally have their 
indices down. 



2 Relativistic wave equations 
2.1 Klein- Gordon equation 

The Klein-Gordon equation was originally thought up by Schrodinger, who 
wanted a relativistic wave equation describing the electron, even before he 
settled on what we now call the Schrodinger equation, 

Klein and Gordon pubhshed first, so they got the credit. We will see later 
why Schrodinger dropped this equation. One of the ways of understanding 
the Schrodinger equation is to recall that in quantum mechanics, physical 
quantities are represented by operators: 

d 

E ^ ih—, p ^ -ihV, (2.2) 

(where the double-headed arrow symbol means "is represented by" ) . Hence 
the Schrodinger equation (2.1) represents the non-relativistic equation for 
the kinetic energy E = p'^/2m. 

As we know, the relativistic relation between energy and momentum is 
£■2 = -I- m^, which seems to suggest that a relativistic version of the 
Schrodinger equation ought to be 

-^0 = -VV + mV, (2.3) 

or, in a manifestly Lorentz covariant form, 

(a' + m^)(j> = 0. (2.4) 

(It is traditional to use (f) rather than ip in this context) . 

Solutions to this equation with definite energy and momentum are 



^g-iBt+p.x ^2.5) 
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with E — ± — -\/(p^ + m^), and A an arbitrary complex constant. 

In fact if is a solution to the Klein- Gordon equation it cannot be 
interpreted as a wave function as Schrodinger discovered. A wave func- 
tion is a probability amplitude, whose modulus squared is the probability 
of finding the particle at a particular position (or with a particular mo- 
mentum). For example, in non-relativistic quantum mechanics, the prob- 
ability density p = which is associated with a probability current 
j = —i{^p*W'ip — 'V'ip*^p)/2m, in the sense that together they make up a 
probability conservation equation 

P + V.j = 0, (2.6) 

an equation which one can check by differentiating p with respect to time 
and using the Schrodinger equation. 

The KG equation also has a conserved density and a 3- vector current, 
which are 

P = —irOfP - dt4>*4>). j = -^(0*V0 - V0V). (2.7) 

However, this density cannot be interpreted as a probability density, as it is 
not positive definite. This was the reason that Schrodinger chose the non- 
relativistic form for his equation. This was reasonable for his purposes, but 
the Klein-Gordon equation makes a comeback later, when we shall sec that p 
can be interpreted as a charge density, which is allowed to take both positive 
and negative values. 



2.2 Maxwell's equations in covariant form 

This section recaps some important results, and introduces the formulation 
of the theory in an explicitly Lorentz covariant manner. 

Firstly, we recall Maxwell's equations in free space, writing them down 
in natural units, in which the permittivity and permeability of free space //q 
and eo are both unity: 

Homogeneous Inhomogeneous 

^ V.B = ^ V.E = p (2 8) 

— B+VAE=0 -— E+VAB=j 

ot at 

The homogeneous and inhomogeneous (i.e. having a source term on the right 
hand side) equations have differing status. The homogeneous equations imply 
the existence of potentials (p and A from which the physically measurable 
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quantities E and B can be calculated. These potentials are specified only 
up to a gauge transformation, — > — A and A — > A + VA, where A is an 
arbitrary function of space and time. The inhomogeneous equations imply 
that the source terms must obey a current conservation equation, /?+ V'j = 0. 
To summarise: 



Homogeneous Inhomogeneous 

potentials 0, A current conservation 

B = VAA 
E = -A - V(/) 

^ 0-A 
A ^ A + VA 



(2.9) 



All these quantities can be assembled into explicitly Lorentz covariant ob- 
jects. The gauge potentials belong together in a 4-vector potential — 
(0, A), while the charge density p and the current density j can be put to- 
gether into a 4-vector current density = (p, j). Recall that putting quanti- 
ties together into 4- vectors is not just a matter of notation: it means that the 
quantities transform just like the space-time coordinates under a Lorentz 
transformation. 

The electric and magnetic fields E and B also belong together in a Lorentz 
covariant object, as they are mixed up by Lorentz transformations (an ob- 
server moving through a magnetic field also sees an electric field). However, 
this object cannot be a 4-vector as there are a total of 6 components of the 
electric and magnetic fields, when they are taken together. In fact, the object 
is an antisymmetric tensor, the field strength tensor F^^'^, which is defined in 
terms of the gauge potential 4-vector A^: 

F^"" = d'^A'' - d^'Ai'. (2.10) 

The Lorentz transformation law for F'^^ follows from those of d'^ and A^^: 

pt^y _^ p'l^^ ^ A'^pA^^F^'^, (2.11) 

from which we explicitly see that F^^ is a covariant tensor of rank 2. Let us 
examine the components where u — Q and ranges over spatial indices i: 



d'A^ - ^A' = d'A^ -A'^ -di(t) -A'^ E\ (2.12) 



where we have used the expression for the electric field in terms of the gauge 
potentials in equation (2.9). The magnetic field is contained in the entries 
where both // and v take spatial values: 

pij ^ Qi^j _ Qjj^i ^ _g.j^j + Q.j^i ^ -Siji^B'', (2.13) 
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where we have introduced the Levi-Civita symbol €ijk- The Levi-Civita sym- 
bol is defined by 



+1 Hi j k cychc, 

-1 iii^ j ^ k anticychc, (2-14) 
otherwise. 



Written in matrix form, the field strength tensor is: 

/ -E^ -E'^ -E^ \ 

The reverse relations may be written 

E' = B' = -^Q,■fcF^■^ (2.16) 

which can be verified with the help of the identity 

^ijk^klm — hl^jm — Sim^jl- (2-17) 

When expressed in terms of 4-vectors and tensors, electromagnetism looks 
very simple and beautiful. For example, the covariant expression of the 
current conservation equation is simply 

= 0. (2.18) 

Maxwell's equations become: 

Homogeneous Inhomogeneous 



(2.19) 



The first of these expressions looks as if it contains many more equations 
than the 4 of the original homogeneous Maxwell equations. However, the 
antisymmetry of the field strength tensor (F^*^ = —P'^p-^ means that the 
expression is trivial if any of the two indices are equal. Thus all the indices 
A, n and u must take different values, and the number of ways of choosing 
three different numbers from a set of four is "^Ca, which is equal to 4, precisely 
the number of equations we started with. 

We can use the four-dimensional Levi-Civita tensor to re-express the ho- 
mogeneous equations more compactly. This tensor has four indices, and is 
defined by 



+1 if 7^ z/ 7^ p 7^ (T, symmetric, 

—1 ii ji ^ u ^ p ^ a, antisymmetric, (2.20) 

otherwise. 
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A symmetric permutation of the indices in one in which an even number of 
indices are exchanged, while an antisymmetric permutation is one for which 
an odd number of indices are exchanged. Thus, for example, 60123 = ^1032 = 

+ I1 ^1023 = ^0132 = — Ij but 60012 = 0. 

There is also a version with the indices raised: 

gMvpa ^ ^t,a^ul3^frr^aS^^^_^^^ ^2.21) 

whose symmetric permutations take the value —1 and antisymmetric ones 
+1. 

Using this totally antisymmetric (i.e. antisymmetric on the exchange of 
any two indices) tensor, the homogeneous Maxwell equations become 

e^^upad'F'^'' = 0. (2.22) 

One of the problems uses this expression to give you practice with the sum- 
mation convention as applied to these tensors. 

2.2.1 Gauges 

Recall that the gauge potential is specified only up to an arbitrary func- 
tion of space and time A (see equation 2.9). The 4- vector version of a gauge 
transformation is written 

A^'^A^'- d^'A, (2.23) 

under which one can check that F^"^ is unchanged. Indeed, using (2.10) one 
finds 

Ff"" d''{A''-d''A)-d%A''-d^A) = F^'' -{d^'d'' -d''d'')A ^ F'"' . (2.24) 

The last step follows because we can take the partial derivatives in any order. 

In order to use the gauge potential to solve the field equations, one should 
eliminate this freedom to make gauge transformations, otherwise one can get 
misled into thinking that the solutions one obtains are all physically distinct. 
There are many ways to do this, but some of the most common are as follows. 

Coulomb gauge This is also called radiation gauge, and is defined by 

V-A = 0, (2.25) 
which implies that the first of the inhomogeneous equations in (2.8) becomes 

= -p. (2.26) 
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This gauge is therefore very useful in solving electrostatics problems, and we 
shall use it later on when studying the Casimir effect. However, it is not so 
often used in relativistic applications, as the gauge condition does not respect 
Lorentz invariancc. 

The second of the inhomogeneous Maxwell equations becomes, on using 
the gauge condition, 

-^(-A- V(/>)-V^A=j, (2.27) 

which we may write as 

(^-V^)A=jT, (2.28) 

where jr = j + V0. One can show, using the equation of current conservation 
and Eq. (2.26), that 

V-jT = 0. (2.29) 

This of course is necessary for the consistency of Eq. (2.28), because if one 
takes the divergence of the left hand side one gets zero as well. 

In free space, that is, in the absence of charges and currents, we may take 
= 0, and the equation for the vector potential becomes 

- V^)A = 0. (2.30) 

One of the important discoveries of the last century was that this equation has 
plane wave solutions, which carry energy and momentum: electromagnetic 
waves. Such a solution may be written 

A(t,x) =ae-^'^*+^'**^ (2.31) 

where a is a constant 3- vector, and u) = ±|k|. 

The gauge condition (2.25) imposes a condition on the vector a, for 

V- A = ik-a e-^'^*+^^"^ = 0, (2.32) 

which implies 

k-a = 0, (2.33) 

Hence the gauge potential A is orthogonal to the wave vector k, which re- 
moves one of the three apparent degrees of freedom of the field. Thus the 
field oscillates in directions orthogonal to the direction of propagation: we 
say that the waves are transverse. 
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Lorentz gauge This is very often the gauge to choose if one is interested 
in wave propagation in electromagnetism. It is defined by 

di.A'' = 0, (2.34) 

and is manifestly Lorentz invariant. In this gauge the equation of motion for 
the gauge field is simplified, for 

^^F^'•' = - d^d'^A^' = = f. (2.35) 

However, an irritating feature of this gauge is that it does not quite specify 
A^^ fully. One can still make a gauge transformation A^ ^ A^ — d^h. which 
satisfies the Lorentz gauge condition (2.34), as long as the function A satisfies 
= 0. Such functions are called harmonic. In classical field theory this is 
not too much of a problem, but in setting up the quantum theory of gauge 
fields care must be taken. 

In free space, the field equation is (2.35) 

□A'^ = - V')^'' = 0, (2.36) 

which looks very much like four copies of the Klein-Gordon equation for fields 
with zero mass. Plane wave solutions may be written 

A''(t,x) = a^e* (2.37) 

where we recall that k-x — — k-x, and — ±|k|. Again, is a constant 
4- vector. 

The Lorentz gauge condition (2.34) are 

d^A^ = -ik^a^e-'^^'''^ = 0, (2.38) 

which implies that 

A; • a = 0. (2.39) 

Once again the field is orthogonal to the wave vector A;'*, but this time in the 
4- vector sense. 

Once the gauge condition is taken into account there are apparently three 
degrees of freedom, one more than in the Coulomb gauge. However, one of 
those is unphysical as a result of the remaining freedom to make harmonic 
gauge transformations. 
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2.3 Lorentz gauge solutions 

The general solution may be constructed from a superposition of plane wave 
solutions ^ 

A^'{x) = / ^ (a''(k)e-''=-^ + a*^(k)e''=-") , (2.40) 

where the amplitudes a'^(k) must all satisfy k ■ a(k) = 0. 

It would be more convenient to be able to choose the amplitudes of the 
components of the gauge field independently, and to this end we introduce a 
set of four basis 4- vectors, or polarisation vectors, e^'^(k), with ^4 = 0, 1, 2, 3, 
for each wavcvcctor k. They must satisfy a completeness condition in order 
for them to be considered as basis vectors, and it is also convenient to make 
them orthonormal, in a 4-vector sense. 

eA*'^(k)eB"(k)77^, = Vab (Orthonormality), (2.41) 
eA*''(k)eB'^(k)77^^ = rj^'' (Completeness). (2.42) 



A good choice is 



eo = (i;0), 

ei = (0, n X k)/|n x k|, 

62 = (0,k X (n X k))/|k X (n X k)|, 

63 = (0;k), 



(2.43) 



where n = (0, 0, 1) is a unit 3- vector. We refer to these polarisation vec- 
tors as timelike (^4 = 0), transverse {A — 1,2), and longitudinal {A — 3) 
respectively. 

Using this basis we can write = a^CA^, and find that the gauge condi- 
tion implies k ■ a = k^a^ — |k|a'^ = 0, or 

a°(k) =a='(k). (2.44) 

Thus the general solution is now more easily written down as 

A^{x) = / ^ (a^(k)eA'^(k)e-''=-^ + a^(k)e^^(k)e''=-^) . (2.45) 

There now appears to be, for each k, three independent degrees of freedom, 
down from the original four. However, one of the amplitudes, a^, does not 
affect the physical value of the electric and magnetic fields. This is a result of 
the remaining freedom in the Lorentz gauge to make gauge transformations 
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with a harmonic function A(x), which may be decomposed into plane waves, 
according to 

A(x) = / ^ (A'^(k)e-''=-^ + A'^(k)e'*^-^) . (2.46) 
Hence under such a gauge transformation, 

a^(k) a'^{k) = a'^(k) - iA;''A(k). (2.47) 
By multiplying with we see that this is equivalent to shifting a^(k), 

a^{k) ^ a'^(k) = a^{k) - i\k\X(k). (2.48) 

Thus although the Lorentz gauge does not entirely remove the freedom to 
make gauge transformations, it does confine the ambiguity in the solution to 
one and only one of the sets of constants in the solution. There are only two 
physical amplitudes to choose, and a^. 

Let us end this section by calculating the Hamiltonian in the Lorentz 
gauge. 

Hi^g = / d^x (vr^i^ - £Lg) 

= -^J d^xiA^.A^ + VA^-VA^'). (2.49) 

We already notice from this expression that the Hamiltonian is not positive 
definite, as the timelike component of the gauge field ^4° appears with a 
negative sign. It turns out that we have already cured this problem in the 
classical theory. To see this, let us substitute the plane wave expansion of 
the field operator, obtaining 

HLg = - f ^u;^VABa*^{k)a^{k). (2.50) 

We saw earlier that the Lorentz gauge condition Eq. 2.34 implies a° = a^, 
and so we find that 




(2.51) 



which is clearly positive definite, and depends only on the physical tranversely 
polarised components. 
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3 Lagrangian formulation of field theory 



3.1 Lagrangian and Hamiltonian mechanics 

The starting point for classical mechanics is Newton's Second Law of Motion, 
which states that the rate of change of momentum of a particle is proportional 
to the applied force. For a system of P particles, we may write 

Pa = Fa, (3.1) 

where the index A runs from 1 to P. The momentum = m^v^, where 
YA — y^A- If the mass is constant in time, we obtain the formulation m^x^ = 
Fa. 

If the force is conservative, which means that the work done by the force 
around a closed path is zero, then the force can be written as the gradient 
of a potential energy V: 

Fa = -VaT^(xi,...,xp). (3.2) 

There may also be constraints on the motion: for example, the particles 
may be forced to move on a sphere. This reduces the number of degrees 
of freedom of the system. A very common form of constraint is one which 
relates some or all of the coordinates through a function (which may also 
depend explicitly on time): 

/(xi,...,xp,i) = 0. (3.3) 

These are known as holonomic constraints. For example, the first particle 
may be constrained to move on a sphere of radius a, or x^ = a^. Thus we 
lose one degree of freedon, the radial coordinate of the particle. If there are k 
such constraints, then the effective number of degrees of freedom is reduced 
by A;: P particles moving in 3 dimensions under k constraints have 3P — k 
degrees of freedom. Thus not all of the coordinates arc needed to describe 
the motion of the system, and it is convenient to introduce the notion of 
generalised coordinates Qa which result from the solution of the constraint 
equations. In the system of particles there will he N — 3P — k oi them, and 
one can express the original coordinates as functions 

Xi = xi(gi, . . . ,gjv,t) 

: ; (3.4) 

Xp = xp(gi, . . . ,gAr,t) 

Generalised coordinates are also useful in systems without constraints: for 
example, orbits around central potentials, where it is convenient to use spher- 
ical polar coordinates. 
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3.1.1 Lagrangian Mechanics 

There is a later and more sophisticated formulation of classical mechanics 
due to Lagrange. We define the kinetic and potential energies T and U in 
terms of the generahsed coordinates of a system as 

U ^U{qaAa) (3.5) 

(we have allowed for the possibility of a velocity-dependent potential). We 
then introduce a new function called the Lagrangian, defined as 

L^T-U. (3.6) 

The equations of motion of the system are then 

These equations are known as the Euler- Lagrange equations. 

For example, suppose there is no velocity dependence in the potential. 
Then 

dL dU dL , , 

^ = ^ = maQa: (3.8) 

SO that 

^(m„g„) = - — . (3.9) 

Thus we return to Newton's Second Law. The quantities dL/dqa can be iden- 
tified as a momentum associated with the coordinate qa - if the qa are simply 
Cartesian coordinates, then they are the ordinary mechanical momenta, but 

in general we refer to them as conjugate or canonical momenta, which are 
defined as the derivative of the Lagrangian with respect to the rate of change 
of the generalised coordinates: 

dL 

Pa = w^. 3.10 
dqa 

The Euler- Lagrange equations can be derived from a variational principle, 
known as Hamilton's Principle or the Principle of Least Action (which is a 
misnomer). The action for a system between times ti and t2 is defined as 



r dtL{qa,qa,t). (3.11) 
Ju 
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The Principle of Least Action states that the action is an extremum for the 
path of the motion. 

To see this, let us consider a small variation in the path, vanishing at ti 
and t2, sends qa(t) — > g^(t) = Qait) + Sqait) (see Figure 3.1). The action also 
changes, I I' = I + 61 , where to first order in 6qa{t), 

61 = dt6L{qa,qa,t) = dt i-Q^Sqa + ^'^^-j " (3-12) 



Integrating by parts we find 



61 



dt (— - — (—]] 6q + 
, \dqa dt [dqa J I 



w-oqa 
dqa 



t2 



(3.13) 



The second term vanishes, as we are keeping the end points of the path fixed, 
and so if the variation of the action 61 vanishes we find that 



dL d_ fdL \ _ ^ 
dqa dt \dqaj 

which are precisely the Euler-Lagrange equations (3.7). 



(3.14) 




q(tj) 

^ q 



Figure 3.1: A path q{t) for a system with one coordinate, and a small 
variation 6q{t) to that path. Classical paths are ones for which the 
action (see Equation 3.11) is stationary under such small variations. 
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3.1.2 Hamiltonian Mechanics 



There is yet another formulation of classical mechanics due to Hamilton. 
It may seem excessive to have so many ways of doing the same thing, but 
each formulation has its own advantages: for example, the Lagrangian and 
Hamiltonian formulations are very powerful in bringing out conservation laws 
in dynamical systems, and relating them to symmetries. 

Recall that associated with each coordinate Qa there is also a conjugate 
momentum Pa, defined by Eq. (3.10). We can therefore express the Euler- 
Lagrange equations as Pa — dL/dqa- We now define an important function 
called the Hamiltonian from the Lagrangian: 

H{Pa,qa,t) =J2Pa<ia-L{qa,qa,t). (3.15) 
a 

The idea is to replace qa by pa, having solved equation (3.10). The first of 
Hamilton's equations follows from the Euler- Lagrange equations (3.7), after 
noting that dL/dqa = —dH/dqa- It is 

dH 

Pa = 3.16 
dqa 

The second follows straightforwardly from the definition of the Hamiltonian 
(3.15) upon partial differentiation with respect to Pa- 

= (3.17) 

dpa 

If the Lagrangian does not depend explicitly on time, the Hamiltonian is a 
constant of the motion, as the following calculation shows. Recall that H is 
a function of ga, Pa, and possibly t only. Hence 

Using Hamilton's equations (3.16,3.17) we see that 

dH_dH__dL 

dt ~ dt ~ df ^^-^^^ 
Thus if dL/dt — 0, the Hamiltonian is constant. One can also show that 

H^T + U, (3.20) 

where T and U are the kinetic and potential energies introduced in (3.5). 
Thus the Hamiltonian can be identified with the total energy of the system. 
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As a trivial example, consider a set of N free particles with mass m. The 
Lagrangian is 

L^Wmql (3.21) 

^ a 

from which we can easily derive the Euler-Lagrange equations 

0-^(mg„) = 0, (3.22) 

or that the acceleration of each particle is zero. The conjugate momentum 
is Pa = mqg, and hence the Hamiltonian becomes 

1 

H^-Y^. (3.23) 
This is of course the familiar non-relativistic expression for the kinetic energy. 



3.2 Lagrangian mechanics for the real scalar field 

A field is an object which takes a value (e.g. a real number) at every point 
in space-time. We are probably most familiar with the electric and magnetic 
vector fields E(t,x) and B(t, x), which we have seen can be unified into the 
electromagnetic field strength tensor F^''{x). The Klein-Gordon field 0(t, x), 
which takes complex values, is another example. 

A very powerful way of classifying fields is to organise them according to 
their properties under Lorentz transformations. The Klein-Gordon field is 
an example of a scalar field, which transforms as a scalar (a pure number): 

<l>{x)^<l>'{x')^cl>{x), (3.24) 

where x'^ — A^^x'^. 

In field theory, the dynamical coordinates, the analogues of the Qa of 
classical Lagrangian mechanics, are the values of the fields at every point. 
A field therefore has an uncountably infinite number of degrees of freedom, 
which is the source of many of the difficulties in field theory. Sums over the 
label a are replaced by integrals over space: for example, the Lagrangian can 
be expressed as the integral over space of an object called the Lagrangian 
density C. For a scalar field, the Lagrangian density can be expressed as a 
function of the field, its time derivative, and also its spatial derivatives: 

L = j d^x C{(f), (j), t). (3.25) 

The appearance of spatial derivatives may at first sight be puzzling, but 
arises naturally from terms coupling neighbouring space points, in the limit 
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that the separation goes to zero. The puzzle is perhaps more that higher 
derivatives don't appear: in practice, they seem not be relevant for most 
applications in particle physics. 

t 



A 




X 



Figure 3.2: A space-time diagram (with the z coordinate sup- 
pressed) of the four- dimensional region in which the variation of 
the field 50(a;) is non-zero. 



From a Lagrangian we can derive Euler-Lagrangc equations, by the ap- 
phcation of Hamilton's Principle to the action for a scalar field. 



Ldt^ dt£x /:{(l),(j),V(l),t). 
Jti Jti 



(3.26) 



Consider an arbitrary variation in the field, (f){x) 4>{x) + S(j){x), where 
S4>{x) vanishes at ti and t2, and outside an arbitrary spatial region R with 
boundary B (see Figure 3.2). This variation causes a small change in the 
action: 

I ^ I' ^1 + 61. 



Taylor expanding the Lagrangian density, we find 

r*2 



r'~ / dtd-'x 
Jti 



dC . , dC , dC , 
90 9V0 d(b 



where 



dC 



dC 



dC 



dC 



dV(t> \d{d,ct>y d{dy4>y d{d. 



(3.27) 



(3.28) 



(3.29) 
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We can separate out the variation in the action 51, and integrate the last 
term by parts, to obtain 



51= r dt(fx 
Jtl 



dC,^ dC 



dt \d4), 







+ 





t2 



Jtl 



(3.30) 



In this last expression the last term vanishes because of the conditions placed 
on 5(j) at ti and t2- Hence, upon integrating by parts again, we get 



51 



rt2 



/ dtd^x 
Jtl 



d_ (dC 



t2 



+ dt dS 



dC 



ti JB dV(p 
(3.31) 

Again, the last term vanishes, as we supposed that the variation 50 vanished 
on the boundary B. Hence we are left with 



51 



t2 



dtd^x 



dC 
dV(f) 



(3.32) 



The argument proceeds by noting that not only is 6(f) arbitrary, but so is the 
region R. Hence the condition that the action be stationary reduces to 



dC 



dC 



d dC 



<9V0 dt d(j) 



0. 



(3.33) 



This is the Euler-Lagrange equation for a scalar field. 

We can put it into a more obviously relativistic form by noting that 
d^(f) = {dtcji-idicj)), so that the Euler-Lagrange equation becomes 



dC dC 

d^ 'm^ 



or 



dn 



dC 



dC 



d{,d^cj>) 



0, 



dC_ 

d(t> 'd{d,4>) 



= 0. 



(3.34) 
(3.35) 



If there are many fields 0a we can derive an Euler-Lagrange equation for each 
of them, 

50a d{df,(f)a) 

We are now in a position to ask the question: from what Lagrangian 
density does the Klein-Gordon equation follow? The Klein-Gordon field may 
be split into its real and imaginary parts, (f){x) = (Pr{x) + i(f)i{x), each of 
which obey the same equation 



{d^d^ + m^)<i)R{x) = 0, {d^d^ + m^)<Pi{x) 



0. 



(3.37) 
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Consider first the following Lagrangian density for the real part, 



We can immediately derive 



d(t) 



dC 



R 



Hence, 



and so 



dC 



d, 



d{di<t>R) 
dC 

dC 



dC 



^R, 



d{dQ(t)R) 



-m^(PR - d^d^ct^R = 0. 



(3.38) 

(3.39) 

(3.40) 
(3.41) 



^d{df,^R) 

Thus the Lagrangian (3.38) implies the correct equation for the real part 
of the Klein-Gordon field: exactly the same argument may be made for the 
imaginary part of the field, so we can write for the total lagrangian 



£ = £r + £ 



= \d^<\>Rd,<\>R-\m^\^\d^<\>id^h-\m^\ 

= \d^<^*d^<^-\m^*<^. (3.42) 



Lastly in this section, we will write down the Hamiltonian for a single 
real scalar field. Drawing on our experience with finite numbers of degrees 
of freedom, it seems natural to define a conjugate momentum 

= ^3.43) 

Thus for our example (3.38), the conjugate momnetum 7r(a;) = 4>{x). The 
Hamiltonian for the real scalar field is defined by 



E = J £x 7r{x)^{x) - J d^xC{(j), di(j), 0). 
If the field obeys the Klein-Gordon equation, its Hamiltonian is 



H 



I 



(fx 



(3.44) 



(3.45) 



We talk of the three terms in the Hamiltonian as kinetic, gradient and po- 
tential terms respectively. One can allow more general functions of (p in the 
potential term than just a quadratic: 



d^x 



(3.46) 
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In classical field theory, V can be any function of 0, as long as it is bounded 
below, otherwise the energy of the field can become arbitrarily negative. 
Quantum field theory is much more restrictive about the terms it allows: it 
turns out that the theory does not make sense in four space-time dimensions 
unless V"(0) is restricted to be a polynomial of degree no more than four. 



3.3 Lagrangian mechanics for the electromagnetic field 

The electromagnetic gauge potential A/^^x) is an example of a 4-vector field, 
which is a field which transforms hke a covariant 4-vector. Its behaviour 
under Lorentz transformations is 

A^ix) ^ A'^ix') = A^A^ix), (3.47) 

where x"^ — A'^x'^. In this section we will write down a Lagrangian density, 
and show that the resulting Euler-Lagrange equations are just the inhomo- 
geneous Maxwell equations. The homogeneous ones will be shown to be an 
identity, rather than arising from any variational principle. The Lagrangian 
is 

jr = ~F,,F'^''-fA,, (3.48) 

where F^i, = d^A^ — d^A^. The generalised coordinates of this system are 
the four components the gauge field A^. We can therefore derive Maxwell's 
equations by looking at the Euler-Lagrange equations: 

^ dC dC ^ 

Let us first consider the variation with respect to Aj,, which gives 

^ = -J'^ = -J'K = -f- (3.50) 



OA, 'OA, ' " 

The variation with respect to d^Ai, is a little more complex, and it will help 
to define X^,, — df^A^,, so that F^^, — X^,^ — X^.^. Hence 

. 9^ ^ 1 dF^ 

"did^A,) idX,^ 4 '"dX^,' 

— _1(^M'^ _ ^^M) _ 1(^/*'^ _ ^i'^^) 

= -F^"". (3.51) 
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Thus we can substitute (3.50) and (3.51) back into the Euler-Lagrange equa- 
tion (3.49) to obtain 

-d^F'^'^+f^O, (3.52) 

which we recognise as the covariant form of the inhomogeneous Maxwell 
equations. 

We now turn to the homogeneous Maxwell equations, and demonstrate 
that they are an identity, which is sometimes known as the Bianchi identity. 
This can be seen by examining the 4-vector — ^e^,^pad'^ F^"' — 0. Sub- 
stituting the expression for the field strength tensor in terms of the gauge 
potential, we obtain 

= le^.p.d''{d'A'' - d^A^). (3.53) 
By exchanging the labels a and p in the second term we find 

= le,,p.d''d^A^-^e,,,,d''d^A'^, (3.54) 
= e^^pad^d'^A'', (3.55) 

with the last step following by the antisymmetry of the Levi-Civita tensor. 
Now we exchange the labels u and p, from which we can deduce 

e^^p^d'-d'^A'' = e^p^.d^d^'A^ (3.56) 

The antisymmetry of the Levi-Civita tensor then shows that 

e^^p.d^'dPA'^ = -e^^p^d^d^A^. (3.57) 

The only number that is equal to its negative is zero, hence 

\e^.pad''F''" = 0. (3.58) 

Thus we have demonstrated that the Principle of Least Action applied to 
the Lagrangian (3.48), combined with the Bianchi identity (3.58), results in 
Maxwell's equations. 

3.4 Noether's theorem and conservation laws 

It is a profound feature of both classical and quantum dynamics that con- 
servation laws, such as conservation of energy or momentum, follow from 
symmetries of the dynamical equations. In this section we shall show how 
this feature is reahsed in classical field theory. 



33 



A conservation law can be expressed in two equivalent ways. Firstly, and 
most obviously, one can say that there is a quantity Q which is constant in 
time, i.e. 

f = 0. ,3.5.) 

Entirely equivalently, one can also say that there is a 4- vector current j'^{x) 
which satisfies a continuity equation, 

d^f{x) = 0. (3.60) 

If we integrate (3.60) over all space, we find 

J £x ^f{x) + J d^xV ■j = 0. (3.61) 

Defining Q — J d^xj'^{x), we find using Gauss's Law that 

-/clS-j, (3.62) 



dQ 
dt 



where the surface integral is taken over a sphere at spatial infinity. With the 
assumption that all currents die off at infinity, we have 

f^O. (3.63) 

Hence the existence of a 4-vector current satisfying a continuity equation 
d^j^{x) = automatically implies that there is a conserved quantity, which 
is the integral over all space of the time component of the current. 

Now consider a theory with N real fields 4>a{x), with a = 1, . . . , N, and 
an action 

S = J d^xC{(t>a,d(t>a). (3.64) 

Noether's theorem for this theory can be stated as follows. 

For every transformation (f)a{x) = <Pa{^') which leaves the action S in- 
variant, there is a conserved current. 

The proof of this theorem is quite involved, and we first need some defi- 
nitions. 

Let us first define the total variation in the field, 

5(l>a{x)^(l>'^{x')-(l>a{x). (3.65) 

Note that this variation can be split into two pieces, the first due to the fact 
that the functions 0a are changing: 

5(j)aix) = (i)'^{x) - (paix), (3.66) 
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which defines what we mean by the operator 5. The remaining piece can 
be interpreted as the change in the field functions due to the change in 
coordinates, 

Maix) = - ~ (3.67) 

where 5x^ — x'^ — x^. Hence 

5(t>a{x) ~ ~5^a{x) + A(/.;(x) ~ ~5(t>a{x) + A(/.,(x), (3.68) 

where we have dropped at term of second order in the variation, as A0„(a;) = 
/^(t)a{x) + 0{5^5x). 

One can show straightforwardly from the definition of 5 that 

~5{d^(t>a{x)) = d^{5(t>a{x)), (3.69) 

but that 

5{d^^a{x)) = d^{54>a{x)) - d,^a{x)d^Sxr (3.70) 

Hence the partial derivative operator commutes with 5, but not with 5. 
Now consider the effect of the variation in Eq. 3.65 on the action: 

5S^ ( d^x'C'ix')- [ d^xCix), (3.71) 
Jn Jci 

where C'{x') = C{(f)'^{x'), d(f)'^{x')), and Q is an arbitrary region in spacetime. 
By adding and subtracting C{x) to the first integrand, we see 

5S ^ j^dt^x' 5C\x) + j^d^x' C{x) - j^d^x C{x). (3.72) 

We now need to express the integration measure d'^x' in terms of d^x, 
using the Jacobian of the coordinate transformation: 



d^'x' = d^x 



ayx ,...,x ) _^4^|^g^M|^ ^373) 



9(x°, . . .,x^) 

where the components of the matrix M are given by 

M^,^-^c^5i^ + dJx'^. (3.74) 

In matrix notation, we may write M = 1 + 5M, and using the standard 
matrix identity 

IndetM^trlnM, (3.75) 

we see that 

detM- l+tr5M = l + ^^^x'^. (3.76) 
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Substituting into the varation of the action Eq. (3.71), we find 

SS^ d^x {l+df,5x^')5£'{x)+J^ d^x {l+df,5x^)£{x)- d^x£{x). (3.77) 

Dropping terms which are second order in the variation we arrive at 

SS^ d^x 5C{x) + ^ d'^x d^dxi'Cix) (3.78) 

Now we can again spht the variation in the Lagrangian function into a piece 
arising from the change in the function itself, and that due to the change of 
coordinates: 

5C{x) = SC{x) + d^C{x)5x''. (3.79) 

The change in the Lagrangian function arises because of the change in the 
functional form of the fields, so 

dC ~ dC ~ 

5C{x) = —{x)5U^) + ^^(x)<5(9>„(x)). (3.80) 

At this point, the commutativity of S and d/j, stated earlier means that 

+ ^ d*x {d^£{x)6x'' + £{x)d^6x^')^ . (3.81) 
Some simple algebra shows that 

SS^ [ d^x (1^ - 9^— 1^ ) 50,(x)+ / d^xd^ { iz^f-rHa + CixjSx"] . 



(3.82) 

Note that the first term is precisely the variation in the action we found 
when deriving the Euler-Lagrange equation, where we considered variations 
in the fields alone, vanishing on the boundary of f2. Suppose we consider a 
particular field configuration (p{x) satisfying the Euler-Lagrange equation. If 
we then make an infinitesimal transformation on it which leaves the action 
invariant, i.e. which results in no first-order change in the action, we must 
have that 

= 0, (3.83) 

where 

dC ~ 

= ^ J 4>a{x) + (3.84) 
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or substituting back for 5(f){x), 

dC ( dL \ 

= ^7^^'^'^«(^) - T^iTHT^^-^'^-^'^^ - (3.85) 
didM \d{d^(j)a) J 

Thus we see the appearance of an infinitesimal 4- vector, f^, which satisfies 
a continuity equation. 



3.4.1 Conservation of energy-momentum 

Consider a translation on the coordinates, with the fields unchanged: 

Sx'' = e", 5(j)a = 0, (3.86) 

where is a constant 4-vector. This clearly includes both space and time 
translations. Hence we can write 

f = r^e^ (3.87) 

where we have introduced the canonical energy-momentum tensor 

dC 

= -wnr-M-^'') - (3-88) 

d[d^j,(pa) 

Because the translations are arbitrary, translation invariance of the La- 
grangian implies that there must be four continuity relations, 

d^e% = 0. (3.89) 

The conserved "charge" associated with these four continuity relations is the 
total 4-momentum, 

= y d'xe\{x). (3.90) 



4 Canonical quantisation of the scalar field 

We have seen that the Klein-Gordon field (j){x) does not have the interpre- 
tation as a single-particle wave-function, as it does not have a conserved 
probabihty current, and it possesses negative energy states which cannot 
form a "Dirac sea". However, as we shall see in this section, the theory 
of the quantum field operator 0(a;) is in fact a relativisitic many-particle 
quantum theory. The process of developing the theory of a quantum field 
is sometimes misleadingly called "second quantisation", perhaps because of 
some idea that one is quantising a wavef unction, itself a quantum object. It 



37 



should be emphasised at the outset that the scalar field, real or complex, is 
not a wavefunction. 

In this section we go through the quantisation procedure for free fields. 
By "free" we mean non-interacting: once a state has been set up there are no 
transitions to any other states. In practice this means that the Lagrangian 
can be written as a quadratic function of the fields. Free fields may sound 
somewhat irrelevant, but nearly all quantum field theories have to be treated 
as small perturbations away from free field theories, so it is an appropriate 
place to start. There are non-trivial effects in free field theory which come 
from boundary conditions: one of these effects is the Casimir effect in the 
free electromagnetic field. 

4.1 Quantisation of finite systems 

In this section we shall proceed with quantising the real scalar field. Firstly, 
we will go through the quantisation procedure for a system with a finite 
number of degrees of freedom, to show that quantising a field is not unlike 
quantising a more familiar system. Suppose the system has coordinates 
Qa with associated velocities Qa, and Lagrangian 

1 ^ 

L^-Y.mql-U{qa) (4.1) 

^ a=l 

(we are supposing for simplicity's sake that all the masses are the same). 
Each coordinate is also associated with a conjugate momentum p^, defined 
as the differential of the Lagrangian with respect to the velocity qa, which in 
this case has the value mqa- The Hamiltonian, defined as | QaPa — L, is 
then 

N 2 

The procedure of canonical quantisation replaces the classical conjugate vari- 
ables qa, Pa by operators obeying the following commutation relations 

[Qa,Pb] = iSab, [qa, Qb] = ^ \Pa,Pb]- (4.3) 

The Hamiltonian is also replaced by an operator H, which generates the 
time evolution of the system. We are probably used to thinking about time 
evolution as the change in time of the state of the system, governed by the 

Schrddinger equation 

ti\m) = H\m)- (4.4) 
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However, there is another way of thinking about the time evolution of a 

quantum mechanical system known as the Heisenberg picture (to distinguish 
it from the Schrddinger picture above). 

In the Heisenberg picture, the states \iP)h that arc time-independent, and 
the time evolution of the system is transferred to the operators Onit), which 
satisfy the Heisenberg equations 

ijOH{t)^[OH{t),H]. (4.5) 

The states and operators in these two pictures can be related by a unitary 
transformation e'^* (a unitary operator is one whose inverse is its Hermitian 
conjugate). Thus, if we suppose that the two sets of operators and states 
coincide at i = 0, they are related at all times by 

Onit) = e'^^'Ose-'"'' 
\i;)H = e^"'\m)s. (4.6) 

This ensures that the matrix elements are the same. To see this, let us 
consider the matrix elements between any two states \iP)h and \'4'')h in the 
Heisenberg picture, and apply these transformations: 

H{i^\OH{t)mH = s{m\e-''''e^'''Ose-^V\i^'{t))s 

= s{m\Os\i^'it))s. (4.7) 

This demonstrates that the Heisenberg picture matrix elements are the same 
as those in the Schrodinger one, which means there is no difference in the 
physical predictions one makes. 



4.2 The real scalar field 

We begin with the Lagrangian density 

1 



£ = -5^09^0 --mV, (4.8) 



from which we can define the conjugate momentum 

Ax) - (4.9) 
d(p[x) 

The Hamiltonian density is then found from 

n = n(x)^(x) - >C = ^TT^ + + \m^\ (4.10) 
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In free quantum field theory we usually work in the Heisenberg picture, 
where the fields carry the time dependence, although it is possible to use 
the Schrodinger picture. The field theory is quantised by constructing the 
equal time canonical commutation relations 

[0(i,x),7r(i,x')] = iM^(x-x'), (4.11) 
[0(i,x),(/)(i,x')] = 0, (4.12) 
[7r(i,x),7r(i,x')] = 0, (4.13) 

where 5^(x — x') = d{x^ — x'^)S{x'^ — x''^)S{x^ — x'^). One can see the analogy 
between the field commutation relations and (4.3) if one recalls that the space 
coordinate x is like the label a in the system with a finite number of degrees 
of freedom. The 5-function acts in the same way for a continuous label x as 
the Kroneckcr 6 does for a discrete label. 

So in quantum field theory, the fields become operators, obeying the 
Heisenberg equations 



1 



1 



0(x) = - (/)(x),// , n(x)^- n(x),H , (4.14) 

where H is the quantum Hamiltonian, constructed by replacing the classical 
fields by operators in the original: 

^-ld'x[ln^ + l{V4>r + l^^'^'). (4.15) 
It is straightforward to show for a free field that 

^{x) = Tr{x), ^x) = (V^ - m^)(i){x). (4.16) 

Hence the field operator obeys the same equation of motion as the classical 
field. It is convenient to expand the field operator in terms of eigenfunctions 
of the operator (— + m^), 

(-V' + m2)e^'**" = c^^e^"^*". (4.17) 

Writing 

<t>{x) = y (/k(x)a(k) + /^(a;)a*(k)) . (4.18) 

with /k(a;) = A]^{t)e^^'^, we see that /k(a;) obeys the field equation (9^ + 
m2)/k(a;) = if 

^k(t) = A^e-''^"* (4.19) 
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The normalisation factor jVk is essentially arbitrary, although almost ev- 
eryone uses one of two conventions outlined below. The properties of the 
operators a(k) and a*(k) can be deduced from the commutation relations of 
(j) and TT. We therefore need expressions for a(k) and a*(k) in terms of and 
TT, where 

7r(x) = '^{x) = j (-^u^,N^,a{\<i)e-"'■^ + ia;kjVk*a*(k)e'*^-^) . (4.20) 

Firstly, we take the Fourier transform of the field operator and its conjugate 
momentum: 

^{k',t) = J £x4){x)e-'^'''', TT{k',t) = J d^xTT{x)e-'^'''', (4.21) 



and use the relation 



J d^xe^^'^-'^')-" = {27rf6%k - k') (4.22) 



to obtain 



4){k',t) = Afka{k')e-'''^'' + Af^a*{-k')e''^^'' (4.23) 
7r(k',t) = -iu;k'A/'k'a(k')e-^'^'''* + ia;k'AAk*a(k')a*(-k')e*"'''*. (4.24) 



Hence 



2uJkAf^^d{k') = J d^x{i7r{x) +uji,4{x))e''''"'. (4.25) 
By complex conjugation one obtains 

2a;kA/'4,a*(k') = J d^x{-m{x) + u;k4{x))e-''''-^ . (4.26) 

One can then derive the commutation relation 

[a(k), a*(k')] = {27:fS%k - k')/2a;k|A4r, (4.27) 

This looks very similar to the commutation relations for simple harmonic 
oscillators. We recall that the harmonic oscillator of unit mass is described 
by the Hamiltonian 

H=Ip' + \u?x\ (4.28) 

where lo is the angular frequency of oscillation. Upon quantisation, the oper- 
ators a = {ip + ijjx)! \f%ji and a* = (— ip-|-a;x)/-s/2a; satisfy the commutation 
relations 

[a,a*] = l. (4.29) 
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This correspondence suggests that one choice for the normahsation factor: 

in which case we find for the complete set of commutation relations 

[a(k),a*(k')] = {27rf6\k-k'), (4.31) 
[a(k),a(k')] = 0, (4.32) 
[a*(k),a*(k')] = 0. (4.33) 

The right hand side of the first equation, with its 5-fucntion, is the closest we 
can get to unity when we are dealing with functions of a continuous variable 
like a(k). However, we shall adopt a convention which keeps the Fourier 
expansions of the field operators Lorentz invariant, 

= ^, (4.34) 

which means that the first commutation relation becomes 

[a(k), a*(k')] = 2iUk{2TTfS%k - k'). (4.35) 

It is also quite common to work in a finite volume V, which means that 
instead of integrals over k we have sums over allowed wavevectors k = 
2'K{ni,n2,n3)/L, which go over to integrals in the infinite volume limit: 

'E-ZtS^- (4-36) 



J (27r) 

All wavefunctions are multiplied by a factor 1 /y/V '■ hence the plane wave 
expansion of the field operator is written 

<^(^) = E {a{k)e-''- + a*{k)e''-) . (4.37) 

The non-trivial equal time canonical commutation relation then becomes 

[a(k),a*(k')]=(5k,k', (4.38) 

where 5k,k' = Sn^n^Sn^n'Jnan'^- 

To summarise, we have discovered that the field theory behaves very 
much like infinitely many harmonic oscillators, labelled by their wavenumber 
k, and with a different frequency tUk = ^/{k'^ + rri^). We will see that the 
occupation number n of each oscillator actually corresponds to the number 
of particles with that particular momentum in the state considered. The 
raising and lowering operators a*(k) and a(k) are also called creation and 
annihilation operators, because they change the number of particles with 
3-momentum k in a state. 
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4.3 States of the scalar field; zero point energy 

In this section we will drop the hat notation for operators, unless there could 
be ambiguity, as all fields will be assumed to be operators unless otherwise 
specified. We saw in the previous section that the scalar field could be 
decomposed into a sum over operators obeying oscillator-like commutation 
relations, one pair for each wavcvcctor k. It will therefore be no surprise to 
see that the states of the quantum field theory are just towers of oscillator 
states, one for each k. The ground state is defined as the state annihilated 
by all the annihilation operators: 

a(k)|0) = 0, Vk. (4.39) 

We can use the creation operators a*(k) to construct excited states: for 
example, the first excited state in mode k may be written 

|k) = a*(k)|0). (4.40) 

The normahsation of these states follows from deciding that the vacuum state 
has unit norm, that is (0|0) = 1. Thus 

(k'|k) = (0|a(k')a*(k)|0) (4.41) 
= (0|a*(k)a(k')|0) + (0|0)(27r)35(k-k'). (4.42) 

The first term on the right hand side vanishes, and so we have that the first 
excited states satisfy the normalisation condition 

(k'|k) = (27r)^5(k-k'). (4.43) 

The corresponding completeness relation is 

d^k 1 



(27r)3 2uk 



|k)(k| = l. (4.44) 



The states |k) are interpreted in quantum field theory as single-particle 
states, an interpretation discussed further below. We could continue ap- 
plying creation operators a*(ki), a*(k2), ... to make states 

|k, ki, k2, ...) = .. . a*(k2)a*(ki)a*(k)|0). (4.45) 

The general state could have multiple applications of creation operators with 
the same momenta, in which case it is conventional to divide by a symmetry 
factor, to obtain 

(n* "l^k 

IW}) = n4=|0), (4.46) 
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where rik is the number of times the creation operator with label k is applied. 
The Hilbert space spanned by the set of all states \nk) is called Fock space. 
The basis states are interpreted as multiparticle states, which demonstrates 
one of the special features of field theory: it is a quantum mechanical theory 
of many particles. Ordinary non-relativistic quantum mechanics describes 
only one particle at a time. 

There is an operator which allows us to count the number of particles in 
a given state. Let us first introduce the set of operators n(k) by 

n(k) = a*(k)a(k). (4.47) 

It is straightforward to see that they return zero when acting on the ground 

state, 

n(k)|0) = 0, (4.48) 

from the definition of the vacuum state (4.39). When acting on 1-particle 
states, we find 

n(k)|k') = a*(k)a(k)a*(k')|0) (4.49) 
= a*(k) (a*(kXk) + 2a;k(27r)3(^(k - k')) |0), (4.50) 

where we have used the commutation relations (4.33) to change the order of 
a(k) and a*(k'). Thus, 

n(k) |k') = 2cuk(27r)35(k - k') |k) . (4.51) 

The operators n(k) are used to construct the number operator 

which tells us the total number of particles in a particular state Thus, as 
expected, A^|0) = 0, meaning that there are no particles in the ground state. 
In the 1-particle state |k'), we find 

^totlk') = / ^^^(k)lk') = / ^(27r)^5(k-k')|k) = |k'). (4.53) 

The state |k') is indeed an eigenstate of the number operator, with eigenvalue 
1 - there is 1 particle in state. When apphed to the general state (4.46), it 
can be shown that 

N\{n^}) - I (03^k|K}). (4.54) 
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The energy of the states can be computed from the quantum energy operator 

E^J Sxe\{x) = / d'x i^-e + \{V^f + , (4.55) 

which is immediately seen to be identical to the Hamiltonian. Since the sim- 
ple states we have seen so far have been constructed from creation operators 

a*(k) acting on the vacuum, it is more convenient to re-express the Hamilto- 
nian operator H in terms of creation and annihilation operators. Recalling 
the expansions of (p and n in terms of a(k) and a*(k) in equations (4.18) and 
(4.20), we find 

^ = /|75^Y(«'(''Mk) + «(k)a-(k)). (4.56) 

If we try to change the order of the operators, we run into a problem, for 

^ ^ / (S^2^T (2«*(k)«(k) + (27r)='5(0)a;k) . (4.57) 

The second term in the brackets has the formally infinite term (27r)^5(0). 
The source of this infinity is the fact that we are working in infinite spatial 
volume. Let us revert to a finite cubic region with sides of length L, in which 
case we can write 



(2nfS(0)^ lim lim / d^xe'""'^^-^^ ^ lim / d^x^V, (4.58) 

L^ook'->-k J_L/2 L-^ooJ-L/2 



where V is a. formal volume factor, and the limit of infinite volume should 
strictly be taken at the end of the calculation of any physical quantity. Thus 
we can interpret the quantity multiplied by V as an energy density. However, 
even this energy density is infinite, in general. Suppose we try to find the 
energy density of the ground state po, for which 

^Vo|0) = j ^c.ka*(k)a(k)|0) + -y a;k(27r)^^c)-(O)|0). 

(4.59) 

Hence, 



This integral is divergent. If we integrate the momenta over the range < 
|k| < A, with the upper cut-off in the integration range A ^ /x, we find that 
the dominant behaviour in the integral is 

= (2^ £ df^f^'if'" + t'")' = (i^^'(l + 0(/^VA'))- (4.61) 
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As we take A to infinity, the ground state energy diverges as the fourth 
power of the cut-off in the momentum - we say that the integral is quariically 
divergent. Thus the ground state energy density of a field theory appears to 
be infinite. This is the first of many infinities in quantum field theory, and 
the way of dealing with this, and all others, is to renormalise, that is, to 
subtract off unobservable quantities. 

As the theory stands, we have no way of measuring the absolute value 
of the the energy of the ground state: all we see is differences between that 
state and others. Thus we measure all energies relative to the ground state 
by subtracting the term which gives the infinite value for the ground state 
energy, | / d^ku!\i5{0), and we define a new, renormalised Hamiltonian 

\ t f d h \ 

H,en^H--j d'ku^SiO) = J ^_a;ka*(k)a(k). (4.62) 

This has the sensible property that Hren\0) — - the ground state energy is 
zero. You might like to check that the energy of the state |k') is indeed u-k': 
a simple way to start is to note that 

f d'k 



'ren 



/d K 
^a;kn(k), (4.63) 



and use Eq. (4.51). 

This conclusion is not sound one we start thinking about gravity. To see 
this, consider the Priedmann equation in cosmology 

= ^p, (4.64) 

where H — d/a is the Hubble expansion rate, a is the scale factor of the 
Universe, p is the total energy density, and G is Newton's constant. This 
equation seems to tell is that we can measure the energy density of the vac- 
uum by measuring the expansion rate of the Universe, and subtracting off the 
contributions from known sources to which we have more direct experimental 
access. In practice, measuring the vacuum energy density (or equivalently a 
cosmological constant) difficult, but through a combination of supernova and 
Cosmic Microwave Background observations, the value expressed in natural 
units is 

PA - 10"^^ GeV^ (4.65) 

Furthermore, it is now generally recognised that the quantum theory of grav- 
ity (if it exists at all) is not a field theory, and that field theories in general 
are only an acceptable description of fundamental particles and forces up 
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to a momentum of around the Planck mass Mp = 1/\/G, and therefore it 
is makes no sense to let the upper limit of integration go beyond Mp in 
Eq. (4.61). There is an analogy here with the theory of phonons in a solid 
lattice: there arc no modes with momenta above 1/a, where a is the spac- 
ing between atoms, but at low momenta we have a very good description of 
phonons as the quantised oscillation in a massless field whose waves moves at 
the sound speed. Hence, in the absence of other interactions besides gravity, 
the vacuum energy of a real scalar field seems to be 

^o^(4^^p. (4-66) 

which is a factor of about 10^^^ larger the the currently accepted value. 

We can recover a more realistic result by starting from a theory with a 
potential 

V{(|>)^Vo + \^x^(t>^ (4.67) 

with Vq — —po + pp^. However, it is quite incredible that nature would have 
arranged such a precisely set value for Vq: po a-nd pa must differ by a precisely 
defined number which is a factor of 10^^^^ smaller than the individual terms. 
This is an example of a "fine-tuning" problem. 

Returning to Minkowski space, and ignoring gravity again, we can sub- 
tract off ground state energies automatically by a procedure known as normal 
ordering. When we place an operator in normal order, we write it so that all 
the creation operators are to the left, and all the annihilation operators are 
to the right, so for example 

:a(k)a*(k'): = a*(kXk). (4.68) 

Hence, 

:H: = //,en. (4.69) 
4.4 Particle interpretation of states 

We shall finish this section by returning to the justification for interpreting 
the states of the field as multiparticle states. Let us first compute the energy 
of a one-particle state: 

i/,en|k) = / a;,,n(k')|k) = u;k|k). (4.70) 

Hence the energy of the state we have decided is a one-particle state is tUk, 
which we recall is equal to A/(k^ + m'^), the energy of a relativistic particle of 
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mass m and momentum k. Thus in order to strengthen our case we should 
check that the momentum of the state |k) is indeed k. 
The momentum operator in field theory is 

P = /d^.^WV^W = /p-^— kf,(k). (4.71) 

There is in fact no need to normal order the momentum operator, as the 
momentum of the vacuum is zero as a result of its translational invariance. 
Thus we may evaluate the momentum of the state |k) through 

Thus the state |k) is an eigenstate of the momentum operator, with eigen- 
value k, properties we expect from a one-particle state. 

The state is also rotationally invariant: a property following from the 
rotational invariance of the scalar field, but should be explicitly checked by 
constructing the angular momentum operators. Hence the state is a spinless 
particle. Lastly, multiparticlc states automatically have the correct Bose 
symmetry. Consider a state of two particles with momenta ki and k2. Then 

|kik2) = a*(ki)a*(k2)|0)a*(k2)a*(ki)|0) = |k2ki), (4.73) 

which shows that the quantum state is symmetric under the exchange of 
particles. This property stems from the canonical commutation relations of 
the field operators. 



5 Canonical quantisation of the electromag- 
netic field 

In order to quantise the electromagnetic field we shall follow exactly the same 
procedure as we did for the real scalar field. We shall find the dynamical 
coordinates, identify the conjugate momentum, and impose canonical com- 
mutation relations between them. However, we should be aware that there 
are extra complications not present for the real scalar field. Firstly, there is 
gauge invariance: recall that the transformation A^j, ^ — d^A does not 
change any physical quantity associated with the field, which means we are 
in danger of encoutering unphysical degrees of freedom when we use A^ as 
the dynamical coordinates. Secondly, the conjugate momentum to A^ looks 
like it should be 

dC 



TT' 



didoA^) 



F^^". (5.1) 
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But this means that does not have a canonical momentum, as tt'^ vanishes 
identically. 

Our approach to dealing with these problems is to choose a gauge from 
the outset, trying to leave only physical degrees of freedom in the Lagrangian. 
We shall carry out this procedure in the Coulomb gauge and in the Lorentz 
gauge. 

5.1 Coulomb gauge quantisation 

We recall that the Coulomb gauge is defined by diA^ — 0. In this gauge, Aq 
is entirely determined at any time by the charge distribution, through the 
equation 

-V^^ = p, (5.2) 

and it is therefore not a true dynamical coordinate. In free space, that is, in 
the absence of charges or currents, we are entitled to set it to zero. 

We recall the action of the electromagnetic field, in the absence of an 
external 4-current: 



S = Jd*x£ = - J d^x ^Ff.^F^", 



where F^,^ = — dvA^. In the Coulomb gauge, this action is equivalent 
to 

5cg = j d^'xCcg = / d'^x (^^A'A' - ^djA'djA^^ . (5.3) 
Thus the momentum conjugate to A'^ is 

dAi ^ ^ 

and the Hamiltonian density is 

^cg = TTiA' - £cg = \a'A' + \djA'd^A\ (5.5) 
At this point it would seem natural to impose a commutation relation ^ 

[A^(x, t), 7r,(x', t)] = iSijSiic - x'), (5.6) 

based on our experience with scalar fields. This is incorrect, as differentiating 
both sides with respect to quickly shows. The gauge condition diA' = 
implies that the left hand side vanishes. However, 

[diA'{x, t), 7rj(x', t)] = idj5{^ - x') ^ 0. (5.7) 



^Note that as this is not a covariant equation, there is no need for covariant and 
contravariant indices to balance on both sides of the equation. 
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The derivative of a 5-f unction is not zero. So we need a more general function 

Aij(x, x') for the right hand side of the CCR, which satisfies 9, A.;, (x, x') = 
(and also 9^ Aij(x, x') = 0, as one can easily check by substituting for iij in 
the commutation relations). We also know that the dependence on x and x' 
must be in the combination x — x' through translational invariance. 
It is easier to find the function Ajj(x — x') in Fourier space: we write 

A,,(x - x') = / ^A,,(k)e^>^-(---'). (5.8) 

Then 

dAiji^ (0^'^^.-(k)e''^-(---') = 0, (5.9) 
from which we infer that 

k'Kij{}i) = 0. (5.10) 
The solution to this equation is 

Chi Kj \ 

- -^j m, (5.11) 

where -F(k) is some function of k. The correct function is just F = 1: this 
keeps the function Ajj(x — x') as close as possible to the 5-function. Thus 
the canonical commutation relations in the Coulomb gauge are 

[^^(x, t), F^(x', t)] = -^ / (0 - ^) g^'^-^^-^'^- (5.12) 
This we can write 

[A^(x, t), E^i^', t)] = -i (s.j - ^) (5(x - x'), (5.13) 

where is the inverse of = didi, and can be thought of as being 
defined in terms of its Fourier transform. Given any function /(x) with 
Fourier transform /(k), V^^/(x) is that function whose Fourier transform is 
k^^/(k). One can now verify that the Coulomb gauge commutation relations 
5.13 are now consistent with the gauge condition diA^ — 0. 

As we saw in Section 2.2.1, the free space Coulomb gauge equations of 
motion for the electromagnetic field are 

{d',-df)A\t,^) = 0, 
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which have the general solution 

where fejO,* = 0. Thus the three components of the vector amplitude a*(k) are 
not independent: they satisfy a constraint arising from the Coulomb gauge 
condition. Rather than worry about this constraint all the time, it is much 
more convenient to have two freely chosen functions of k, so we instead write 
the general solution 

r d'^h 1 

- Y.J (2^72^ [h^Ax)^-'-' + -WU--'-') ■ (5-14) 

The e^;^) are called polarisation vectors, and have the properties 

fc'e(A) = 0, e\^/(^^,^ = 5xx'. (5.15) 

These conditions mean that, along with fc* = A';V|k|, the two polarisation 
vectors form an orthonormal basis for 3-vectors. The completeness property 
therefore follows: 

^(1)^(1) + ^(2)^(2) + ^'^■^ — ^ij^ (5.16) 

or EAe(A)e(A) = k - ■ 

Upon quantisation, the independent amplitudes a(A)(k) become opera- 
tors, and by substituting the expansions for A'(i, k) and E^{t,yi') into the 
canonical commutation relations, and using the properties of the polarisation 
vectors, one discovers the commutation relations for a(A) (k) and a^^^ (k) : 

[a(,)(k), a^,,)(k')] = 5M27rf5(k - k'). (5.17) 

Thus each polarisation behaves like an independent scalar field. 
As before, we can define a vacuum state |0) satisfying 

a(;,)(k)|0) = 0, Vk,A. (5.18) 

Other states are derived from the vacuum by acting with raising operators. 
For example, there is a set of "first excited states" 

|k,A) = a^,)(k)|0). (5.19) 

Note that there is an extra label, A, which indicates the polarisation of the 
electromagnetic field in the state. The quantised states of the electromagnetic 
field are of course the particles known as photons. 
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Last in this section, we write down the Hamiltonian for the electromag- 
netic field in terms of the ladder operators in its normal-ordered form: 

/t d h 
"^'"^'^^ ^ J (2^^^'^«(A)(k)«w(k). (5.20) 

This is just two copies of the renormalised scalar field Hamiltonian (4.62), 
one for each polarisation. 

5.2 Canonical quantisation in the Lorentz gauge 

In order to quantise this theory in the canonical manner, we must first iden- 
tify the canonical momentum conjugate to the degrees of freedom, impose 
canonical commuation relations, identify the ladder operators and find their 
commutation relations, and construct a Fock space of states. We should also 
calculate the Hamiltonian, and verify that it is positive definite. 

As one might imagine, gauge invariance poses problems for the quantum 
theory as well. Firstly, we have the problem mentioned in Section 2.2.1, that 
the canonical momentum vanishes. This turns out to be cured by fixing the 
gauge, as we can verify in the Lorentz gauge: 

= ^ = -A,. (5.21) 

Hence the equal time canonical commutation relations can be guessed to be 
those for four independent real fields 

[A^(t,x),7r.(t,x')] =<5-^(x-x'), (5.22) 
[A>'{t, x), A'^it, x')] = 0, [n^it, x), 7r,(t, x')] = 0. (5.23) 

By substituting the plane wave expansion for the field operator A^[x), one 
can find (after some algebra) that 

[a^(k),a*^(k')] = -^7^^2cuk(5'(k-k'), (5.24) 
[a^(k), a^(k')] = 0, [a*^(k), a*^(k')] = 0. (5.25) 

Hence we have four pairs of ladder operators, one pair for each polarisation. 
We define the vacuum state |0) by 

a"^(k)|0) = Vk,A. (5.26) 

and we can construct excited states by acting with the raising operators 
a*^(k). For example, 

|k,/l) = a*^(k)|0). (5.27) 



52 



How do we impose the gauge constraint in the quantum theory? Firstly, 
it is clear that demanding d - A = is too strong a condition on the operator 
A^, as it is inconsistent with the equal time canonical commutation relations 
Eq. 5.22. A weaker condition is to ask that matrix elements of the gauge 
condition vanish. This we cannot do for every possible state, but we only 
accept as physical states for which this is true: i.e. 

• ^10) = (5.28) 

for any two physical states |0) and Any equivalent condition is to 

demand that the positive frequency part of the gauge condition gives zero 
when acting on a physical state: 

d-A^+^\(f))^0, (5.29) 

where ^ 

d ■ = - 1 '^zk ■ ae- ''^-\ (5.30) 

Hence physical states must satisfy 

(a°(k) -a3(k))|(/.) = 0, (5.31) 

which is the equivalent of Eq. (2.44) in the quantum theory. 

Finally, we shall check that the Hamiltonian is positive definite for phys- 
ical states. We must deal with the infinity in the vacuum energy by normal 
ordering, and study the matrix elements of the operator 

:i^Lg: = - / |^a;k77ABa*^(k)a^(k). (5.32) 

The quantum version of the gauge condition Eq. (5.31) implies that for phys- 
ical states 10) and |0'), 

(0'|:i/Lg:|0) = E /f^a;k(0la*''(k)a^(k)|0). (5.33) 

In particular, the expectation value of the energy is clearly positive definite 
for physical states. 



6 Fermions 

Dirac knew of the problem with the probability interpretation of the Klein- 
Gordon equation, and traced it to the fact that it was second order in time 
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derivatives. He therefore set out to find a relativistic wave equation with 
only one time derivative. The requirement that the form of the equation be 
unchanged under Lorentz transformations, which mix up d/dt and V, means 
that the equation must be first order in spatial derivatives as well. Hence 
Dirac proposed a relativistic free particle wave equation 

d 

i— ■0 = —icx-'Vih + dmih. (6.1) 
at 

There clearly has to be something rather special about the objects a, (3 and 
if) in order that Lorentz covariance be preserved. In fact, a and (3 are 4x4 
Hermitean matrices, and ip is an object with four components called a spinor. 

Now, we want a wave equation with plane wave solutions which satisfy 
the relativistic energy-momentum relation E'^ — "p^ + m^. To see how this 
property emerges from the Dirac equation, we act on both sides of (6.1) with 
{idt — ioL • V), to obtain 

( d \ ( d \ ( d \ 

U— — ia-Vj ( i— + ia • V 1 = m ( i— — • V 1 Z^'^. (6.2) 

At this point, it is more convenient to switch to index notation for the 3- 
vectors a and V, noting that a • V = a^di. We expand out the brackets on 
the left hand side of (6.2), and add and subtract —i(3a''di on the right hand 
side. The result is 



(6.3) 



It is now useful to define the anticommutator, represented by curly brackets: 

{A, B}^AB + BA. (6.4) 

(Sometimes you will also see the anticommutator written as [A, 5]+). We 
rewrite the left hand side as (— 9| + \a^aWidj + d^djdi)'^, which we are 
entitled to do by renaming the indices on a and 9, to obtain 

+ a^]did)j ij = m^p'ip + a^jdiij;. (6.5) 

In order to reproduce the relativistic relation between energy and momentum, 
we must end up with an equation like the Klein-Gordon equation. If a* and 
P satisfy 

{a\a^}^S'^, {P,a'}^0, = 1, (6.6) 
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where 1 is the 4x4 identity matrix, t/j satisfies the equation 

Hence, each of the four components of ip satisfy the Klein-Gordon equation. 

We can make the Dirac equation more exphcitly relativistic by defining 
four new 4x4 matrices: 

7° ^(3, y = (3a\ (6.8) 
If we multiply both sides of the Dirac equation (6.1) by P, we obtain 

(i-f^d,, - ml)'ijj = 0. (6.9) 
The conditions (6.6) on a* and P are neatly unified into the matrix equation 

{7/^,7-} = 2?7'^"1. (6.10) 

6.1 Plane wave solutions 

We know that plane wave solutions exist, but what are they? We therefore 
look for solutions of the form 

V^(x) = -ue-'P-^, (6.11) 

with u a constant 4-component spinor. Although it appears from (6.7) that 
u is arbitrary, this is not in fact true, as we shall see. 

Substituting our plane wave ansatz (6.11) into the Dirac equation (6.9), 
we find 

(7% - m)ue~'P-'' = 0. (6.12) 

At this point we shall introduce another piece of notation: we shall use the 
"slash" symbol / to denote contraction of a 4- vector with the Dirac gamma 
matrices: 

0^7^^=^. (6.13) 

Thus u satisfies the equation 

{^-m)u^O. (6.14) 
Let us premultiply this equation by + m), and use an important identity, 

^^=a-61, (6.15) 
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to show that 

{•j^ -\- rn){'^ — rn)u — {p^ — m^)u — (6.16) 

This brings us back to the relation E'^ = + m^, although by a rather more 
rapid route. It is conventional to choose = E always to be positive, in 
which case we must also include solutions proportional to e"*"*^'^ separately. 
We write them 

^(x) = -ye**"^, (6.17) 

which satisfy 

{■^ + m)v^O. (6.18) 

Solutions with time dependence e~*^* (e*^*) are called positive (negative) 
energy solutions. 

We can use the first equality in Eq. (6.16) to show that the following 
spinors are solutions to the Dirac equations for u and v: 

u±{p) = N{^ + m) j v±{p) = -N{^-m) j 

where is a normalisation factor, and we have defined x+ ^^^^ X- to span 
the space of 2-component spinors (which means they are spinor spanners): 

X+=(J), ^-=(?)- (6.19) 

Any solution to the Dirac equation can be written as a superposition of the 
four four-component basis spinors u± and v±. 

To proceed further we need an explicit representation for the gamma 
matrices, and we will use the so-called "standard representation" given in 
Appendix A.l. It will turn out to be convenient to choose the normalisation 
factor to be 

N = l/y/{E + m). (6.20) 

The basis for the 4-component positive energy solutions can therefore be 
written 



(^ + ^)^±V (6.21) 



y/E + m V p • <yx± J ' 

and the negative energy basis is 

y/ETT^ \ iE + m)x± J ^ ' 

There arc a couple of remarks that are worth making about the solutions 
(6.21) and (6.22). Firstly, when the 3-momentum is small, that is when 
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|p| <S m, two of the components are much smaller than the other two. For the 
positive energy solutions, it is the lower two components that are negligible, 
while the opposite is true for the negative energy solutions. Secondly, the 
appearance of two degrees of freedom for each set of solutions (that is, x±) 
demands some explanation. One can verify that x± eigenstates of the 
Pauli matrix cr^, with eigenvalues ±1: 

(t'x± = ±X±- (6.23) 

This degeneracy is actually a result of the fact that the particles described 
by the Dirac equation (such as the electron) have an extra property, spin, 
which is a type of angular momentum intrinsic to the particle itself, unlike 
ordinary angular momentum which particles gain by virtue of rotation. Spin 
is quantised in half- integer units, with the spin of the electron being s — ^. 
As with other types of quantised angular momentum, the value of the spin 
projected along any axis (such as the z axis) takes the values s, s — 1, . . . , —s. 
Thus Dirac-type particles have only two spin states, which we conventionally 
call "up" and "down". 

With the given normalisation convention, the spinors have the properties 

u\uB = 2E6ab, v\vb = 2E6ab (6.24) 
We define the adjoint spinors u and v by 

UA = ^7°, VA = ^;l7°- (6.25) 
The adjoint spinors have the further properties 

uaUb — ^tuSab, vaVb = —2mdAB- (6.26) 

Note that this is a different normalisation convention than that used by 
Greiner and Reinhardt. 



6.2 The interpretation of negative energy states 

In the last section we called half of the solutions "negative energy" solutions 
without proper explanation. Of course, the explanation is that they appear 
to represent quantum states with negative energy, for if one acts with the 
quantum mechanical energy operator, one obtains a negative number: 

i^ive'P-'') = -E{ve'P-''), (6.27) 

with E = |(p^ + m^)2 1. This raises an ugly possibility: that an electron with 
positive energy Ei could decay into one with negative energy —E2, radiating 
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a photon with energy E1 + E2. There indeed seems to be nothing to stop the 
energy of the electron from tumbling down towards minus infinity, releasing 
an infinite amount of energy in the process. This is clearly nonsense. Dirac's 
way of circumventing this problem was to suppose that all the negative en- 
ergy states are filled already, and the Pauli exclusion principle prevents a 
positive energy electron from radiating a photon and occupying a filled nega- 
tive energy state. One can still imagine exciting one of these negative energy 
electrons into a postive energy state, whereupon it becomes a real electron, 
leaving behind a hole. A hole is an absence of a negative electric charge, 
which has positive charge. Thus Dirac realised that his theory predicted 
the existence of a positively charged spin | particle with exactly the same 
mass as the electron. Initially, Dirac identified this particle with the proton, 
hoping that somehow Coulomb interactions would provide a mass difference, 
but after others (including Weyl, Oppenheimcr, and Tamm) showed that this 
couldn't be right, he eventually had to abandon this conservative position 
and predict a new particle now know as the positron. Its discovery in 1932 
by Carl Anderson in cosmic rays settled the issue. 

However, this picture of a negative energy "sea" of electrons leaves much 
to be desired. Bosons, particles with integer spin such as the photon (spin 
1) and the pion (spin 0) do not obey the Pauli Exclusion Principle and 
yet they still have negative energy states. Thus there can be no Dirac sea 
picture and we are left where we were before. It is only if we abandon simple 
wave mechanics and turn to quantum field theory that we find a satisfactory 
resolution of this problem. 

6.3 Quantising the spinor field 

The Lagrangian density for a massive spinor field is 



where wc recall that in the standard representation of the Dirac gamma 
matrices the adjoint spinor = ■ip'^'-y^. Rather like the complex scalar field, 
one can obtain the field equation (which is the Dirac equation) by treating ip 
and -0 as independent quantities, and demand that the action be stationary 
with respect to arbitrary variations in either. The variation of the action 
with respect to 'ip gives 



jC — iip ^ijj — mijjijj, 



(6.28) 




(6.29) 



(the Lagrangian is not a function of d/^ip in this formulation) . Hence 




(6.30) 
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from which we derive the Dirac equation i ^ijj — rmjj — 0. Similarly, we can 
vary with respect to ip and obtain 



SS = j (fx (# ^5ip - mip6ip). (6.31) 
Integrating by parts we find 

5S = j d'^x {-id^i)^ - mi))5^l: + j dS/^iipYS^, (6.32) 

where the last term is an integral over the space-time surface at spatial 
infinity (|x| — > oo) with end-caps at |t| — > oo. As usual, we suppose that the 
variations die away at infinity so that we can drop the surface term, so we 
recover the equation for the adjoint spinor id^ip'f'^ + mtjj = 0. 
The momentum conjugate to ip is given by 

7r = -j= zV^t", (6.33) 
dip 

which shows that in the standard reresentation ip and iip^ are canonically 
conjugate variables. Thus we can find the Hamiltonian density: 

n^TTip-jC^ -ii^-f^dki^ + mi^i^. (6.34) 

The quantisation of the spinor field follows a familiar pattern. We first of 
all suppose that ^^{x) is an operator satisfying some commutation relations, 
acting on some space of quantum states. We specify the equal time com- 
mutation relations, and then try to find the possible states. The problem 
with the spinor field is that it is not obvious from the outset what commu- 
tation relations to impose on the field operator, and the obvious relation 
[^{t, 'K),iip^{t, x')] = iS{x — x') is not in fact correct. 

To see why this is so, let us compute the zero-point energy of the field. 
Firstly, we expand the field in terms of its operator- valued Fourier coefficients: 

^(^) = E / -^3^ {cA{p)uA{p)e-^'- + dl{p)vA{p)e''-) , (6.35) 



A=±- 



where = + m^, and ma(p) and Va{p) are 4-component spinors which 
satisfy 

- m)uA{p) ^ 0, {^ + m)vA{p)^0. (6.36) 
Similarly, the conjugate spinor has the expansion 

V'(^) = E / T^TJ^ (c'Aip)uA{p)e+''-^ + dA{p)vA{p)e-'''-'') . (6.37) 



A=±- 
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We can now compute the Hamiltonian, which is the spatial integral of the 
Hamiltonian density or 

j £x {-iipYdi'ip + mi;^) (6.38) 

The end result is 

^ = E / (0^2^ ^ (cIi(p)ca(p) - ciA(p)4(p)) ■ (6.39) 

The vacuum state is traditionally defined as the state annihilated by all 
annihilation operators ca{p) and ^^(p): 

ca(p)|0) = 0, ciA(p)|0) = 0. (6.40) 

We expect an infinite zero point energy again, removed by normal ordering. 
However, H is not positive definite, and nor is :H: unless we define the 
commutators of the operators such that 

:cZ^(p)4(p): = -4(p)d^(p). (6.41) 

This means that the fermionic creation and annihilation operators must an- 
ticommute, and we quantise the spinor field by anticommutation relations 

{rf^(p),4(p')} = 2EdABHp-p'), (6.42) 
{ca(p),4(p')} = 2ESab3{p-p'). (6.43) 

All other anticommutator brackets vanish: 

{ca{p),cb{p')}= ={dA{p),dB{p')} (6.44) 
{ci(p),cUp')}= ={4(p),4(p')} (6.45) 
{ca(p),4(p')}= ={c\{p),dB{p')}. (6.46) 

There are four kinds of 1-particlc states in this theory, created by the four 
operators c^(p) and 4(p)- We write them 

|p,^;0) = c^i(p)|0), |0;p,^) = 4(p)|0). (6.47) 

As you can probably guess from the way the states are built up for the 
complex scalar field, 0^4 (p) creates a particle with momentum p and spin 
state A, while 4(p) creates an antipariicle. 
Let us examine the 2-particle state 

|pi,Ai,p2,A2;0) = 4(pi)4(p2)|0). (6.48) 
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If we interchange the particles, we must interchange the creation operators, 
which gives the state a relative minus sign: 

|p2,A2,pi,Ai;0) = 4(p2)4(pi)|0), (6.49) 

= -4(Pi)4(P2)|0), (6.50) 
= -|pi,Ai,p2,A2;0). (6.51) 

Thus we see something that has to be taken for granted in ordinary quan- 
tum mechanics: the states of the spinor or Dirac field are antisymmetric 
under particle interchange. We call these states fermionic, and the particles 
fermions. Furthermore, suppose we take both momenta and both spin states 
equal to p and A respectively: we then find 

|p, A p, A; 0) = -|p. A, p. A; 0) = 0. (6.52) 

Hence two fermions can never be in the same state. This is precisely the 
Pauli exclusion principle, which is fundamental to the understanding of the 
physics of the atom. On the other hands, bosons, like the scalar particle we 
considered previously, have commuting field operators, and so any number 
of them can be in the same state. 

Armed with the anticommutation relations for the creation and annihila- 
tion operators, and the spinor completeness relations 

^a^b = (7 • P)ab + m5ab, ^a^b = (T " P)ab " m5ab, (6.53) 

A A 

the fields can be shown to satisfy the equal time anticommutation relations 

{Mt,^),Mt,^')} = iM„;,5(x-x'), (6.54) 
{V'a(t,x),V6(i,x')}= ^{Mt,yi),Mt,^')} (6-55) 

where the factor of h has been revived, and labels a and b have been put in 
to keep track of the rows and columns of the spinors and -i/j. 



6.4 Conserved charge 

One can show straightforwardly from the Dirac equation (and its adjoint) 
that the 4- vector quantity 

j^(a;) = i,{xy/''^{x) (6.56) 
satisfies a continuity relation d ■ j = 0. Hence there is a conserved charge 

Q = / d=^x'07V= / d^xV'V- (6.57) 
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By substituting the plane wave expansion for the Dirac field operator, one 
can show that 

^3 = E/ (0^2^ (4(p)ca(p) - 4(p)^a(p)) . (6.58) 

Applying this operator to to the single particle states, one can show that 

Q\p,A;0)^+\p,A;0), Q\0,p,A)^-\0;p,A). (6.59) 

Hence the states created by c^(p) have equal and opposite charge to those 
created by rf^(p), which fits in with their interpretation as particles and 
antiparticles. 

One can also show that the existence of this conserved charge is a conse- 
quence of a symmetry of the Lagrangian, as we should expect from Noether's 
theorem. This symmetry is 

ilj{x) i/j'{x) = e'^i/jix), i/j{x) xIj'{x) = e-'^^i/jix). (6.60) 

7 Scalar fields in an expanding Universe 

In this section we shall look at a quantum field theory in an expanding 
space-time, as is appropriate for early Universe cosmology, and in particular 
for the generation of perturbations in an infiationary phase. We shall follow 
our normal procedure of starting with an action, deriving the Eulcr-Lagrange 
equations, finding the solutions, and promoting the field and its canonical 
conjugate to operators. We will find the fascinating result that in a universe 
expanding exponentially fast the vacuum is no longer quite empty: the field 
fluctuates with an amplitude H, where H is the expansion rate. In the 
inflationary universe scenario is the traces of these fluctuations which we see 
in the Cosmic Microwave Background. 

7.1 Scalar field in curved spacetime 

In order to define a field theory in a general spacetime background wc must 
go beyond Special Relativity and make it consistent with the postulates of 
General Relativity. This means that we must be able to formulate it in such 
as way that the equations of the theory take the same form for any observer, 
not just those in a state of uniform motion. Thus the action must be invariant 
under a general coordinate transformation 

x'' ^x"" ^x'^ix). (7.1) 
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This is a much larger group of transformations than the Lorentz group of 
Special Relativity. Under this transformation, a rank 2 tensor like the metric 
transforms as 

^ g'^,{x') = 9pa{x)^^^^^^ (7.2) 

It is not too hard to show that the following expression for the action of a 

real scalar field is indeed invariant under general coordinate transformations, 
and reduces to the Minkowski space action of the previous sections when 

S = j d'x^g {^-g^^d,(t>d.(t> - ^(0)) , (7-3) 

where g = det (jf^^ . A flat homogeneous and isotropic expanding Universe is 
described by the metric g^^ = diag(l, — a^(t), — a^(t), — a^(t)). By substitut- 
ing this metric in the action (7.3) we quickly see that it is 

S^l dV (i^^ - ^(V0)^ - y (</>)) . (7.4) 

Using the Euler-Lagrange equation we find 

+ 3-0 - + V'{<j)) = 0. (7.5) 

The scale factor a{t) is determined from the Priedmann equation 



a 

where G is Newton's gravitational constant and p is the energy density. Let 
us suppose that the energy density is dominated by the energy density of a 
scalar field, which is 

P=^0' + ^^(V0)^ + m- (7.7) 

This must obey the covariant form of the energy-momentum conservation 
equation 

p + 3H{p + p)=0 (7.8) 

where H is the Hubble parameter H = a/a, p is the average energy density, 
and p is the average of the pressure p, given by 

P=^0^-^^(V0)^-m- (7.9) 

Let us look for solutions which consist of small fluctuations around a homo- 
geneous field 0(t, x) = (j){t) + 99(t,x), with \(p\ <C |0|. By substituting this 
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Ansatz into the field equation (7.5), and neglecting the fluctuations, we find 
for the homogeneous part: 



+ 3-0 + 1/'(0) =0. (7.10) 
a 

Note that when the is large, the field is strongly damped, so we look for 
solutions in which we can neglect in comparison to 3H^, so that 

3H}c^-V'{^). (7.11) 

One can show that sufficient conditions for this to be true are that the fol- 
lowing two dimensionless quantities are small: 

— — . n — mfi — — 



where rup — l/SnG is the square of the so-called reduced Planck mass. 

These are the so-called slow roll parameters. In the limit e 0, one finds 
that = 0, and that p = —p = V{(f)). Thus the energy density and pressure 
are constant and dominated by the potential energy density V, and solving 
the Priedmann equation one finds that the scale factor expands exponentially: 

a(xexp{Ht), (7.13) 



with H — ^V/3mp. This rapid expansion is an example of inflation. More 
generally, we define infiation to be a situation where a > 0, and we find that 
e is non-zero. In this case the expansion goes approximately as a power-law 
a oc We can regard exponential expansion as the limit of power-law 

expansion in the limit e ^ 0. 

7.2 Classical scalar field fiuctuations 

We now examine the equation for the fluctuations (/?. By keeping only terms 
flrst order in ip in the full field equation (7.5) we find that the fiuctuations 
obey 

^ + 3-^- \v^^ + V"{4>)^ = 0. (7.14) 
Using the slow- roll parameter rj, 

+ 3-0 - 4tVV + Sryi/V - 0, (7.15) 
a 
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where H — a/ a is the Hubble parameter. One definition of an infiating 
universe is that a > (in these coordinates) . One can show that an inflating 
universe whose energy density is dominated by the potential of a scalar field 
V should generally have r] 1. 

Let us now look for solutions to Eq. (7.15). First it is convenient to 
introduce conformal time through the relation adr — dt, or 



f dt'la{t') 



(7.16) 



If a oc t^, we have that r oc ^/(l — p), and hence a oc r^^*^^ ^\ By a 
suitable choice of origin, one can also show that 

- ^ ' (7.17) 



1 — p aH 



Note that exponential expansion corresponds to the limit p — cxd, in which 
case r = —{aH)~^. Hence —dot < 0, and the universe inflates as r — > 0^. 
The general solution to Eq. 7.15 has the form 

Vit, ^) = / 2 {aMt)e'^- + atme-''- '') (7.18) 
where = -\/(/c^ + SrjH'^), and the functions /k satisfy 

A + 3-/k + ^/k + SrjH^h = (7.19) 



a 



Let us now suppose that the Universe is inflating and we can neglect terms 
of order rj. Then to simply the equation we can write fk{t) — a{t)u-k:{t), and 
obtain the equation for Uk{t) 

-Mk+ Uk = 0. (7.20) 



dr^ V T 

One can easily show by direct substitution that the two independent solutions 
are 

f^^t) = i^lfAe-^'^r ^ _H^t^e-'^r_ (721) 
akr k 

and its complex conjugate /k(^)- 

Fluctuations in the scalar held lead to fluctuations in the energy density 
5 p. Because the energy density is dominated by the potential, the flrst order 
term in the fluctuation is 

5p{t, x) = p(t, x) - p = V'{4>)^{t, x). (7.22) 
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The mean square fluctuations in the energy density are therefore 

5p%^)) = {V'm'{^%^)) = ^p'{^%^)). (7.23) 

We will see that in an inflationary universe, the fluctuations are non-zero, 
even in the vacuum state. 



7.3 Quantum scalar field fluctuations 

Let us now quantise the fluctuation in the fleld (p in the familiar canonical 
manner, assuming exponential expansion, and zero mass. The conjugate 
momentum is identified as 

(7.24) 

giving a Hamiltonian 

n = 7r<P-C = a' + ^(V^)2 + y(^)) , (7.25) 

which we note is proportional to but not identical to the energy density. We 
now promote </? and tt to operators and impose the canonical commutation 
relations 

[ip{t, x), 7r(t, x')] = i5^{x - x') (7.26) 

The wave functions /k and are already normalised so that the creation and 
annihilation operators satisfy the standard relativistic commutation relations 

[ok, Ok'] = 2a;p ^(k - k'), [ok, Ok'] = [al., a^,] = 0. (7.27) 

With these creation and annihilation operators we can define the states in 
the same way as in Minkowski space. The vacuum and 1-particles states |0) 
and |k) are defined 

OklO) = Vk, |k) = al\0) (7.28) 

respectively. 

Let us now consider the mean square fluctuations of the fleld in the vac- 
uum state, 

a'^{0\ip\t,^)\0), (7.29) 

which we have seen determines the fluctuations in the energy density. By 
substituting the expansion of the fleld in terms of its ladder operators, and 
using their commutation relations, one flnds 

.^./^|/,(,)p./f(|)\l + *V). (7.30) 
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Define the power spectrum of the field fluctuations to be 




(7.31) 



we see that 

As the universe inflates, r — ^> from negative values, 
second term for an individual mode. Hence 

^.(fc)-(^)'- (7.33) 

Thus we conclude that in an exponentially expanding Universe, a massless 
scalar field fiuctuatcs on all length scales with an amplitude H/2t:. Defining 
the density fiuctuation on a comoving scale k by 

Sk^- [ £x6p(t,x)e-'^'', (7.34) 
pj 

we recover the classic result 




(7.35) 



Strictly speaking, in pure de Sitter space, the energy density is constant and 
= 0. This result should be understood as being the leading term in an 
expansion in the slow roll parameters rj and e. Further more, in order to 
calculate the effect of these energy density perturbations today, one needs to 
evolve them forward through the end of infiation and past their eventually 
re-entry through the horizon, when once again k = (aH). This requires cos- 
mological perturbation theory which is discussed in Martin Kunz's lectures 
in the current volume. 

Based on our experience with the vacuum energy, we may be concerned 
that the result for the mean square fiuctuations (7.30) contains divergent 
quantities: the first term on the right hand side is logarithmically divergent, 
and the second term quadratically. 

Let us examine the divergence in more detail. We again regularise the 
integrals to make algebraic manipulations meaningful. We will again impose 
an upper cut-off A on the magnitude of the physical momentum variable k/a. 
The cut-off is supposed to arise from a parent theory at higher energies, and 
may even be the Planck mass mp . Because of the logarithmic nature of the 
divergence in the first term, we must also consider a lower cut-off on the 



(7.32) 

and we can neglect the 
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integral. This is actually supplied by the effective mass squared = 3r]H 
in the more physical case rj 0. 
Hence the regularised version is 



a^- i^] log (7-37) 



<r..= ^\Mt)\' = T-,, (7.38) 



Using the relation Eq. 7.17 in the limit p — > oo, 

2 /AX A2 

We note that the second term is independent of H, and therefore we should 
suspect that is also present in Minkowski space. Let us denote the fluc- 
tuations in Minkowski space by ctq, where the wave-functions are fk{t) = 
exp(— ia;kt), 

Hence we see that the second term in Eq. (7.36) is identical to the Minkowski 
space fluctuations. The fact that o"^ and ctq diverge should give us pause for 
thought, but is not fatal to the theory. We should bear in mind that the 
upper cut-off cannot be consistently taken higher than mp, and in any case 
we do not observe the fluctuations in the field at a single point in space: we 
always average over some region, giving a perfectly finite result. 

To discuss the effect of these ffuctuations further would take us too far 
afield: we would end up having to consider field theories with terms higher 
than quadratic in the potential, for which we need to develop the machinery 
of interacting field theory. 



8 Interacting fields 1 

So far we have been dealing with free field theories, which we have been able 
to solve exactly. Thus we have found the equations of motion for the field 
operators, solved them in terms of a sum of operator-valued Fourier modes 
(a(k), 6(k) and their complex conjugates) and found all the possible quantum 
states. The dynamics of a free field theory are therefore rather trivial: one 
decides which state (or superposition of states) the field is in at some initial 
time, and the field then remains in this state for all subsequent times. 

We would like to able to describe more realistic situations in which the 
field changes its state. Field theory was developed to describe the result 
of scattering experiments in which a few particles head towards each other 
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t=-T 



t=+T 




Figure 8.1: A scattering experiment. A set of particles with mo- 
menta ki, k2 come in from a large separation a long time in the 
past, interact, and head out towards infinity again in a different 
state, where they have different momenta. There may also be a 
different number of particles of various species in the final state, 
infinity, interact. 



from a long distance, interact, and the resulting particles (there may be more 
than the incoming set) then head back to infinity again (see Figure 8.1). The 
goal is to calculate the quantum mechanical amplitude for the initial state 
to change in to the final state. From this we can calculate the transition 
probability, which is usually expressed as a cross-section. 

Unfortunately, it is almost always impossible to calculate these probabili- 
ties exactly, and one must use a perturbation expansion in some small number 
which parametrises the strength of the interaction - a coupling constant. For 
example, q, the charge of the complex scalar field of the previous section, 
is such an expansion parameter, which we must assume to be small. Most 
of the technical complexity of quantum field theory is due to the difficulty 
of writing this expansion in a manageable way, and in then in making sense 
of the expressions once they have been written down. The breakthrough 
that made field theory calculations feasible was Feyman's realisation that 
the perturbation expansion could be written down in a graphical way with 
Feynman diagrams. The procedure of making sense of the result, which is 
formally infinite, is renormalisation. 
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8.1 The interaction picture 

In order to develop this expansion we must first recall some results from ordi- 
nary time-dependent perturbation theory in quantum mechanics. There are 
three ways of looking at the time evolution of a state in quantum mechanics, 
which are equivalent in the sense that all operators (which correspond to 
observable quantities) have the same matrix elements. These ways are called 
"pictures" . 

Schrodinger picture This is probably the most familiar picture. The time 

dependence of the system is carried in the states, while the operators 
Os are time- independent. The time-dependent states \ip{t))s obey the 
Schrodinger equation 

i^^m))s-H\ms, (8.1) 

where H is the Hamiltonian. 

Heisenberg picture In this picture, it is the states that are time- 

independent, and the time evolution of the system is transferred to the 
operators Onit), which satisfy the Heisenberg equations 

tjOH{t) = [OH{t),H]. (8.2) 

The states and operators in these two pictures can be related by a 

unitary transformation e*'^* (a unitary operator is one whose inverse 
is its hermitean conjugate). Thus, if we suppose that the two sets of 
operators and states coincide at t = 0, they are related at all times by 

\^)h = ^e'"'m))s. (8.3) 

This ensures that the matrix elements are the same. To see this, let us 
consider the matrix elements between any two states \iP)h and \iP')h 
in the Heisenberg picture, and apply these transformations: 

H{i'\OH{m')H = 5(^(t)|e-'^*e^^*05e-^^*e^^*|V^'(t))5 

= s{m\Os\^'{t))s. (8.4) 

This demonstrates the point. 
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Interaction picture This picture is the one suited for situations in which 
we can't solve the Schrodinger (or Heisenberg) equation exactly for our 
Hamiltonian iJ, but we can split it into two parts, one solvable with 
Hamiltonian Hq, and the other H' which is small compared to Hq. 

H^Ho + H'. (8.5) 

In the interaction picture, the time dependence is divided between the 
operators and the states. The operators evolve according to the Heisen- 
berg equations for the solvable Hamiltonian Hq, 

Oi{t) = e^^°*Ose-^^°*, (8.6) 

which means that the states must be related to the Schrodinger picture 
states by 

i±\^{t))j^e^"'''\m)s. (8.7) 
Thus by differentiating both sides with respect to time, 

(8.8) 

The piece is brackets o the right hand side is just the interaction picture 
representation of the interaction Hamiltonian, which we write H'(t)i. 
Thus 

if^m))i-H'(t):m))i- (8.9) 

Thus the interaction picture states evolve according to the non-trivial 
part of the Hamiltonian only. In the limit if' — > 0, the interaction 
picture reduces to the Hesinberg picture. 

The interaction picture state evolution equation has a formal solution 
which is 

m))j^u{t,to)mto))i, (8.10) 

where U{t,to) is an operator which satisfies the equation 

i^U{t,to)^H'{t)iU{t,to), (8.11) 

with the boundary condition f/(t,to) = 1, the identity operator. We can 
convert this into an integral equation 



U{t, to)^l-t [' dhH'j{t,)U{ti, to), (8.12) 
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Figure 8.2: The simplification of the integration region in Equa- 
tion (8.14). The orginal integration is over the darker shaded 
area, but by using the symmetry of the integrand and careful 
time ordering the region can be extended over the entire shaded 
area. 



which we solve by iteration. The starting point is to take U{t,to) — 1. We 
substitute this back into the integral equation, from which we obtain 

U^^\t,to) = l-i f dtiH'j{ti). (8.13) 

Substituting again, we get 

Jto 

= 1 + {-if f dh r dt2H'i{ti)H'j{t2). (8.14) 

Jto Jto 

We can simplify the integration region by exploiting the symmetry of the 
integrand, and by using the time ordered product introduced in Section (8.4). 
Note first that we are integrating over the darker shaded area in Figure (8.2). 
Let us now rewrite the second order term in Equation (8.14) as 

tjl! fdh r dhH'i{t,)H'i{t2) + ^^ fdh r dhH'j{t2)H'j{t,). (8.15) 

Z Jto Jto Z Jto Jto 



72 



All we have done here is to divide the integration into two and interchange 
the integration variables in the second term. In doing so, we can see that 
the second term is effectively an integration over the lighter shaded area. In 
both terms the earlier time appears to the right. Thus we may write this as 

fdh fdhT[H'j{h)H'j{t2)]. (8.16) 
z JtQ J to 

It is not too hard to satsify oneself that the nth term in this iteration involves 
the time ordered product of n copies of the Hamiltonian, with a symmetry 
factor of l/n\. Hence 

U{t,to) = J2^-^ fdh... f dtr,T[H',{ti)...H[{tn)]. (8.17) 
„ n\ Jto Jto 

This can be written more compactly as 
U{t,to)^T 

where the operation of time ordering on an exponential is defined by its 
operation on the individual terms in the Taylor series. 

8.2 The S-matrix and transition amplitudes 

In principle, the matrix C/(t,to), whose formal solution we found in the pre- 
vious section, completely determines the time evolution of the system. In 

scattering processes, we are generally not interested in the intermediate times 
when the particles are interacting, only in how the states in the distant past 
change into other states in the future. We assume that interactions are lo- 
calised in space and time, which amounts to assuming that Hj{t) — > as 
t ±oo. Hence 

\^{t))i ^ |*(±oo)),, (8.19) 
which are states of the free Hamiltonian Hq. We define the S-matrix to be 

lim U(t,to), (8.20) 

so that 

|*(+oo))/ = 5|*(-oo))/. (8.21) 

We call states |^(— oo))/ in states, and |\&(+oo))/ out states. In order to 
reduce the amount of typing, they will be denoted \i) and |/) respectively. 

In a scattering experiment, we prepare the system in a state |^'(t))/ at 
t — > — oo, and we would like to be able to compute the amplitude for a 



exp 



: / dt'H'iit" 

Jtn 



(8.18) 
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quantum mechanical transition to any other state \^{t)) as i — > +00. The 
probabihty amphtude for this transition is 

A^^= hm imt)\^ (t)) I =1 m+oo)\S\^ {-00)) I. (8.22) 

Simphfying the notation, the probabilty amphtude for the transition from 
an in state \i) to an out-state |/) is 

Af, = {f\S\t) ^ Sf,. (8.23) 

That is, the matrix elements of the S-matrix give the transition probability 
amplitudes between initial and final states. 

The S-matrix has a number of very important properties, the first of which 
is that probability must be conserved. That is, if we prepare an in-state with 
unit probability, the sum of all the probabilities of the out-states must also 
be unity. Thus we are assuming that = 1, and that this normalisation 
is preserved throughout the evolution: that is, /(^(-|-oo)|^(-|-oo))/ = 1 also. 
Hence 

{i\S^S\i) = 1. (8.24) 
Let us insert a complete set of out-states: 

l = El/)(/l' (8-25) 
/ 

so that 

j:{t\Sy){f\S\z)^l, (8.26) 
/ 

or 

JlSlfSfi^l. (8.27) 

/ 

Thus the conservation of probability demands that the S-matrix be unitary 
{S^S=1). 

Finally for this section it is important to point out a major problem with 
this approach to dealing with interactin field theories. The interaction picture 
Hamiltonian Hj{t) does not become negligible as t ^ ±cxo, because particles 
are always interacting with themselves. The result of these self-interactions 
is to shift the mass and charge of the particle by an infinite amount, although 
this only happens at second order in the perturbation expansion for S. These 
infinite shifts must be dealt with by the procedure of renormalisation. 
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8.3 Example: X(f)^ 

We will consider a scalar Lagrangian densoty supplemented with a quartic 
interaction term, with a strength parametrised by a constant A: 

C = 1902 - imV' - ^A0^ (8.28) 

Given that the action S = J d'^xC is dimcnsionless in Natural Units (it has 
the same dimensions as h), and that the partial derivative has dimensions 
of inverse length or mass M, we can immediately establish that the scalar 
field has dimension M and that the parameter A is dimensionless. Parameters 
such as A multiplying non-quadratic monomials in field theory lagrangians are 
known as coupling constants, and they measure the strength of interactions 
between particles. 

Prom this Lagrangian density we can easily see that the Hamiltonian 
density is 

7i = Itt^ + + \m'ct>' + (8.29) 

A convenient decomposition into a solvable part and perturbation is to take 

Ho^j d'x Qn' + + lm'<t>') ^ H,^ j d'x ^A^l (8.30) 

Hence the -S-operator in this theory is 



T 



exp (^-i J d^x ^X(l>j{x] 



(8.31) 



with the interaction field operator satisfying 



t^(Pjix)^\(l>iix),Ho\. (8.32) 

Prom now on we will drop the subscript on the interaction picture fields, and 
it will be assumed that all field operators are in the interaction picture unless 
otherwise specified. 

Thus in order to find S-matrix elements we will need to evaluate expres- 
sions containing long products of field operators in time order. We also have 
the option to normal-order the interaction Hamiltonian - the advantage be- 
ing that the vacuum state is left unaffected by the S-operator. Pinding an 
algorithm to evaluate these expressions is the subject of the next few sec- 
tions, which will lead towards a convenient graphical representations known 
as Feynman diagrams (or graphs). 
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8.4 Feynman propagator for scalar field 

A special role is played by the vacuum expectation value of the time-ordered 
product of two fields, which is termed the Feynman propagator: 

iA^ix, - X2) = (O|r[0(xi)0(x2)] |0). (8.33) 

Note that it is translational invariance which tells is that the function Ap{xi,X2) 
can only depend on the difference of its arguments. Prom the definition of 
time ordering, 

iAp{xi - X2) = (0| {(I){x,)cj){x2)e{h - t2) + (t){x2)(t){xi)d{t2 - t,)) |0). (8.34) 

Next we decompose the field operator into positive and negative frequency 
parts: 

0(a;) = 0W(a;) + (j)^~\x), (8.35) 



with 



cl>^^\x) = I ^ay^e-^'-, 0(-)(x) = / ^aie^^'-. (8.36) 



2a;k J 2u}k 

Note that (f)^~\x) — {(f)'^'^\x))'^ . The positive frequency part annihilates the 
vacuum, and so 

iA^ix, - X2) = (O|0W(xi)0(-)(a;2)|O)^(ti - ts) + m'-+\x2)(p'^-\x^)\0)9it2 - h), 

= {o0-'\x,),<i>^~\x2)]\o)e{h-t2) + {Q0+\x2)<i>^~\x,)]\o)e{t2-h). 

(8.37) 



We define a new function 



AM(x,-X2) = -^[0(+^:^l),0^-H^2)]= /^f^k,aj,]e-'^-^+^'='-^ 
^ -if |!^e-*(^i-^^), (8.38) 

and its partner 

A^-\xi-X2) = -i[4>^-\xi),4>^+\x2)] = {A^+\x^-X2))^. (8.39) 

Hence 

iAp^xi - X2) = iA^^\xi - X2)9{h - t2) + iA^-\xi - X2)9{t2 - h). (8.40) 
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8.5 Wick's theorem 

Wick's theorem re-expresses time-ordered products of fields in terms of normal- 
ordered products (whose matrix elements are much simpler to evaluate) and 
Feynman propagators. 

In the following, we will use the notation (pi = (f){xi) for economy. Hence 

cpih - 0S+V5^^ + 0S"V5^^ + 0S^V^"^ + (8.41) 

which differs from the normal-ordered version only in the third term: 

:0i02: = ^S+V^^^ + + 0^"VS+^ + <f>[~Ui~^- (8.42) 

Hence 

0i02-:0i02: = [0^\0^"^], (8.43) 

or 

0102 = :0i02: + iAW(xi - X2). (8.44) 
Armed with Eq. (8.44) we can rewrite the time-ordered product of two fields: 

T[M2] = M20(h-t2) + MAt2-ti), 

= (:(/.i(/.2: + iAW(xi-X2))^(ti-t2) + (1^2). (8.45) 

Given that :4>i4>2'- = •4>24>i'-, we find 

^[0102] = :0i02: + i^F{xi - X2) (8.46) 

that is, the time ordered product of two fields is just the sum of the normal- 
ordered product and the Feynman propagator, or the vacuum expectation 
value of the time-ordered product of the two fields. 

We will use a special notation to denote the vacuum expectation value of 
the time-ordered product, the "tie" or "contraction": 

Hxi)(p{x2) = (O|T[0(xi)0(x2)]|O) (8.47) 

Wick's theorem is the extension of Eq. (8.46) to many operators, and is 
loosely stated by saying: 

The time ordered product of any number of fields is the sum 
of terms containing all possible contractions (including zero) with 
the remaining operators in normal order. 

Note that the contraction of a pair of operators already in normal order 
is zero: 

:0(xi)0(x2): = O. (8.48) 
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8.6 Using Wick's Theorem 

As we know, transition amplitudes follow from S-matrix elements, which 
can be evaluated in the Dyson- Wick expansion. If we want the transition 
amplitude between states l'^) and \phi), the expresstion is 



out(0|-S|V')in=E'^^"H'/''V^) (8-49) 

n=0 

with 

5(")(0,V^) = ^^ / (0|r[:Hi(xi):...:7ii(xJ:]|V^)iA...(iV, (8.50) 
nl Jout 

where the interaction Hamiltonian for our real scalar field theory is given by 

n,{x) = ^X(l>\x). (8.51) 

Hence 

SM{^,^) = t^(^\ f (</,|T[:0^(xi):...:0Vn):]|V'M'a;i...rfV. 
nl y4!y ^out 

(8.52) 

Thus we see that the Dyson- Wick expansion can be viewed as an expansion 
in powers of A. 

We will consider transition amplitudes between multiparticle momentum 
eigenstates |ki, . . . , k„), which obey orthonormailty conditions 

(k; . . . k:„|ki . . . k„) = 6mn2u;y., . . . 2u;^„ Yl - Ki) ■ ■ ■ " KJ^ 

{pi} 

(8.53) 

where {pi} is a permutation of 1 . . .n. That is, the states have zero inner 
product unless the number of particles is the same and the momenta are the 
same. This is easily demonstrated from the definition of the multiparticle 
states in terms of creation operators aj^, acting on the vacuum state |0). 
We will also need the following identities, 

[0(+)(a;),aj,] =e-^- [0(+)(x), ak] = 0, . . 

[0(-na^),«k] = -e*^ [0(-)(a;),4] = O. ^ " ^ 

To be specific, we will consider the S-matrix elements for 2-particle in-states 
and out-states: 

|V')in = |kik2)in = aj^^aj^jO) 

|0)out = |k;k'2)out = 4;a;,,|0). (8.55) 
We now proceed to evaluate S^'"^ for n up to 2. 
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8.6.1 

This is the trivial S-matrix element, as it has no powers of the coupling 
constant and therefore is identical to the free field theory: 

- (k;k'2|kik2) = 2uj^,2uj^, ( 5(ki - k;) 5(k2 - k^) + 5(ki - k'2) 5(k2 - k;)) 

(8.56) 

The particles go "straight through" without interacting. Note that there are 
two terms because the particles are indistinguishable. 

8.6.2 

Now we have the possibility of something interesting happening. Two parti- 
cles may emerge with different momenta, which is the definition of scattering. 

S^'^ = -z^ j out(k;m:0^(a:i):|kik2)in(i'xi (8.57) 

Let us expand the matrix element of the normal-ordered field operator mono- 
mial, dropping the "in" and "out" notation to simplify the expression. In 
scattering amplitudes, this should be unambiguous as in-states are always on 
the right and out-states always on the left. 

(k;k^|:0^(xO:|kik2) = + K^^iW^^f 

+ 6(0S-V(</>i+y + 4(</>i-))^^(0S+)) + (0i-))^] |kik2). 

(8.58) 

Consider now the action of the positive frequency part of the field operator 
on the in-state: 

0W(xi)|kik2) = [0^+H^i),«l4]|o) 

= ([0(+n^i),4j4. + 4j0^+n^i),4j) |o) 

= |k2)e-''^i-"i + |ki)e-*^2-"i. (8.59) 
This can be wriiten in a "tie" notation: 

0W(xi)|kik2) = 4>''+\xi)\kr\^2) + <P^^\xi)\ki\i2) (8.60) 

I I I I 

Thus we learn that when a positive frequency part of a field operator 0^"^^ [x] 
is tied to a particle with momentum k in the in-state, both are removed and 
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replaced by a wave function factor e The operation is clearly repeatable, 
so 

(/.W(a^l)|klk2) = (f>^+\xi)(t>^+\x2)\k,k2)+(l>^^\x2)\kik2) 

I !^ I I I 

= (/)(+)(a;i)0(+)(a;2)|kik2) +0(+)(a;i)0(+)(a^2)|kik2) 
I ' I I I ' 

_ ^^-iki-X2-ik2-xi _j_ ^-iki-xi—ik2-X2^ g^) 

By taking the Hermitian conjugate of these expressions, we learn that we 
can "tie" negative frequency parts to out-states: 

(...k'...|0(-)(a;) = (...|e'^-^ (8.62) 

I I 

Hence 

{kik2\(t>'^-\xi)(t>'^~\x2) = (0| {e^K-^^+iK-^^ + g^'^i'^i+^'^i-^^j . (8.63) 

A very important point to notice is that once all particles in the in-state 
have been tied to postive frequency operators, the vacuum state remains, so 
the application of a further positive frequency operator annihilates the state. 
Thus, in an expression where the field operators have been expanded in terms 
of their positive and negative frequency parts 

• the contribution to the matrix element vanishes if there are more pos- 
itive frequency operators than in-state particles. 

A similar statement obviously exists for the out-state: 

• the contribution to the matrix element vanishes if there are more neg- 
ative frequency operators than out-state particles. 

If there are fewer operators than particles, then once they have all been 
tied off, we will be left with inner products between multiparticle states, 
corresponding again to particles which travel "straight through" without in- 
teracting. These terms will not generally be relevant for scattering. 

Hence we conclude that the only term in Eq. (8.58) which does not vanish 

is 

(k;k'2|:0^(xi):|kik2) = (k'.k'^l [6(0i-))2(0S+))^] |kik2), 

— 4|gj(fei+A;i-A;i-fc2)-a:i_ (8.64) 

Hence 

S^^) ^-iXd{k[ + k'2-ki-k2). (8.65) 
There are two important points to notice here: 
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1. Overall conservation of 4- momentum has appeared in the 5-function, 
as a result of integrating the wave function factors over all space-time. 

2. The 4! factor has cancelled. This happens because there are 4! ways 
to tie the positive and negative frequency operators to external states 
(which was of course the reason that the factor was included in the first 
place) . 



81 



A Dirac matrices 

Dirac matrices (in four space-time dimensions) are 4x4 matrices defined by 
tlie anticommutation relations 

{7^7"} = 2<''l4, (A.l) 

where ■q^'^ is the Minkowski metric, and I4 is the 4x4 identity matrix (which 
is often dropped and left implicit in the equation). 
Properties of the Dirac 7-matrices include: 

(i). (7^)'!" = 7°7'^7°. 
One can also define another, linearly independent, 7-matrix 

f = ijVl'l', (A.2) 
which has the easily derivable properties 

(7^)^ = 1, (7^)t = 7V7°, {7^7^ = 0. (A.3) 

The eigenvector of this matrix are known as chirality eigenvectors. Eigenvec- 
tors with eigenvalue +1 are termed right-handed, with -1 left-handed. One 
can define projectors onto these eigenvectors: 

Pn^li^ + l'), Pl-1{1-i') (A.4) 
A.l Standard representation of Dirac matrices 

Here, (T* arc the Pauli matrices, which are defined below, and and 1 are the 
2x2 zero and identity matrices respectively. 

A.2 Pauli matrices 

The Pauli matrices {i — 1,2,3) are defined to be 

^-(?;). ^^=(roO' ''-(I-.)- 

They are Hermitean, i.e. 

(a^)^ = a\ (A.7) 
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and they square to unity: 

{a^f = (aY = {a^ = 1. (A.8) 

One can verify that 

crV^ = i(7^, (andcydic). (A.9) 
Lastly, any two different Pauh matrices anticommute, or 

(tV^ + (tV^ = 0, (and cychc). (A. 10) 

This set of relations can be neatly expressed using the 3 dimensional Levi- 
Civita symbol e'^^'^: 

[a\a^] = 2ie'^''a\ (A. 11) 
and with the anticommutation relations 

{a\a^} = 25'^. (A. 12) 

One may thus write the product of any two Pauli matrices as 

aV^' = l{a\ a^} + ha\ a^] = 6'^ + ie^^V^ (A.13) 

The Pauli matrices form a representation of the angular momentum al- 
gebra with (spin) angular momentum |. 
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Problem Sheets 



.1 Problem Sheet 1 

1. In the following, is a 4- vector, and = 

(a) Show by explicit different at ion that 

and hence that 
Evaluate d ■ x 

(b) Show that, if 4){x) = Ae~*^-^, with p a constant 4- vector, then 
(9^ + m^)(j){x) — (— + m^)(j){x). 

1 2 

(c) Given the function f{x) = e~^°'^ , where a is a constant scalar, 
and X is a 4- vector coordinate, calculate d^^d^fix) and show that 
d^f{x) vanishes on the curve x'^ = 4/a. 

2. (a) A'^i,, the tensor representing a Lorentz transformation on a four- 

vector, is defined by the property 

Show that if and b" are any two contravariant four-vectors, 
then this property ensures that the scalar product a • 6 is invariant 
under a Lorentz transformation. 

(b) Let (A-^)'*^ be the inverse of A'*^, such that {A-'^^^A'^p = S^. 
Show from the transformation law x ^ x'^ — A^^^.x^ that 



where 9^ = 
(c) Show that 



(A-^)^^ = A/. 



(d) Show that for a two-index tensor M^y that M^i^ and Mi,^ arc not 
necessarily equal (i.e. that the ordering of indices is important). In 
the lectures the Kronecker delta was written 5^, not distinguishing 
between and - why do we not need to worry about index 
ordering in this case? 
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(a) Let ■0(x, t) be a wave function satisfying the Schrodinger equation 

Show that the probabihty density p = {ipl'^ and the probabihty 
current j = —ih{ip*'Vip — ip'Vip*)/2m satisfy a continuity equation 

P+V.j = 0. 

(b) Let 0(a;) be a complex-valued field satisfying the Klein-Gordon 
equation 

+ m^ff) = 0. 
Show that the 4-vector current 

satisfies a covariant continuity equation (9 ■ J = 0. Give a reason 
why can not be a probability (4-)current. 

A real scalar field has Lagrangian density 

where the potential energy density 1^(0) is bounded from below. 

(a) Calculate the canonical energy-momentum tensor Q^^ from this 
Lagrangian density, using the formula derived in the lectures. 

(b) Calculate the Hamiltonian density 7i, and verify that it is equal 
to the energy density 6'"° [x] . 

(c) Verify, using the field equation -|- V\i>) = 0, that the energy- 
momentum tensor is conserved (satisfies 8^9^'^ — 0). 
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B.2 Problem Sheet 2 



1. (a) Write down the Lorentz transformation which takes the 4- vector 

= {ujk, /c*), where ujk = |(k^ + m^)2 1, to a frame moving with 
velocity v in the x direction relative to the original. 

(b) By applying the Lorentz transformation of the previous part, or 
otherwise, show that the momentum integration measure 
is Lorentz invariant, 

2. Consider a real scalar field whose Lagrangian density is C — l^^^^'*^ — 
|m^0^. Find the Hamiltonian operator H, and show that the Heisen- 



berg equations of motion 



7t(x) — - 7t(x), H 



imply that (5^ + m?)^{x) — 0. 
3. The Hamiltonian for a free real scalar field is H — |(7r^ + V0^ + m^0^). 

(a) Show that 

j Sx m'^(t?{x) = 



/ 



(2^ 



(b) Calculate also / d^x (V0(a;))^ and / cPx (j)^{x), and using your 
results, show that 

_ j d^k cuk 



H 



OkOk + OkOk 



(2u;k) 2 
4. Let /k(a;) = e"''''^. 

(a) Show that the functions f-k{x) obey orthogonality conditions 

i j d'x if:ix)doMx) - dof:{x) f^ix)) = 2uj^S'{k - k'), 
i j d^x {fk{x)dofi^>{x) - dofk{x) fk>{x)) = 0. 

(b) Show that, for a real scalar field, 

ai,^i j Sx (^f^{x)n{x) - dof^{x)^{x)) , 
where ^{x) = do(f){x) is the canonical momentum operator. 
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(c) Show that the non-trivial commutation relation for the ladder 
operators is 

5. The number operator for a real scalar field is defined as 

where and Ok are the ladder operators, which obey the commutation 
relation of Problem 4. 

(a) Show that [N , (aj,)"] = n(4)"- 

(b) Hence, or otherwise, show that the state 




is an eigenstate of the number operator with eigenvalue 

where y is a formal volume factor. 
NB The product over the continuous variable k is defined as 




where k = 27rm/L, and m = {m}, m^, m^), with w}, w?, integers. 
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B.3 Problem Sheet 3 

1. In the lectures we defined the functions 

with fco = u;k = ^(k^ + m^). 

(a) Show that A(+)(a;) has the following properties: 

A^+\x') = A(+)(a;), if x"' = AV, 
AW(-,t) = M+\x), 
A^+\x) = 0, if x^<0. 

Here, A^i, is a general Lorentz transform. 

(b) Find A{x - x') = (O|[0(a;), ^{x')]\0) in terms of A^ and A(-). 

2. Consider the one-parameter family of Lagrangians 

(a) Find the field equations for the gauge potential A^^ which result 
from this family. 

(b) Show that they all give the same field equations in the Lorentz 
gauge. What is special about the field equations when C = 1? 

(c) Show that the Lagrangian 

C = --d.A''d''A, 
2 

gives the same field equations as Ci, and find L — L\. Why is it 
that L and L\ give the same field equations? 
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3. (a) Show that the Lagrangian L — —^d^A'^d^A^ gives the Lorentz 

gauge field equations for a free electromagnetic field. Show that 
the associated Hamiltonian density is 

2 ^ 2 ^ 

(b) Given the expansion of the electromagnetic field operator in a 
plane wave basis, 

A,(x) = (a^(k)e^(k)e-*^ + ai(k)et^^(k)e*) , 

with — cok, show that 

:H: = / d^xiH: - - / P^r]^''a^A{k)aB{k). 
J J 2(x^k 

(c) Show that if |0) and arc physical states (in the meaning of the 
Gupta-Bleuler quantisation procedure) then 

{^\:H:\^) = J ^u;k(^|(al(k)ai(k) + 4(k)a2(k))|0). 

4. In the following, (j){x) is a real scalar field operator. 

(a) Show that :4>{xi)4>{x2)'- = :(p{x2)4>{xi): and : 4>{xi)4>{x2) ■ = 0. 

(b) Use Wick's theorem to expand T[0(xi)0(x2)0(x3)0(x4)] in terms 
of contractions and normal ordered products. 

(c) Expand also T[:0(xi)0(x2)::0(x3)0(x4):] 

(d) How many ways are there of taking c contractions of n operators? 
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B.4 Problem Sheet 4 

1. Consider the second order contribution to the S-matrix element for 
2^2 scattering of real scalar particles 



2 



^£2 - J ^'^i^'^2(k;k^|T[:0^(xi)::0^(x2):|kik2) 

(see e.g. Greiner & Reinhardt). Evaluate this expression and show that 
there arc three separate contributions iM^X2 (where 5'£2 — ^-^£2 5^(A;i+ 
k'2-ki- k2) ) 



1 


(p2 


- + ie)[{ki + k2 
1 


-pY 


— m? + ie] 


(p2 


— m? + ie)[{k'i — ki 
1 


-pY 


— + ie] 


(p2 


— + ie)[(/ci — k2 


-pf 


— + ie] 



Write down the Feynman graphs representing these expressions. 

2. Consider a real scalar field theory with potential V{(j)) = |m^0^ + ^A0^. 
A particular scattering process has two particles in the in-state, with 
3-momenta ki, k2, and four particles in the out-state, with 3-momenta 
k'l and k2, kg, and li^. 

The yS-matrix element is written 

54^2 = (k;,k^,k^,k;|5|ki,k2). 

(a) Write down the Dyson- Wick expansion for the S operator, and 
show that the ©(A*^) and O(A^) terms in the corresponding ex- 
pansion of the matrix element vanish. 

(b) Apply Wick's theorem to T[:(f)^{xi)::(f)'^{x2):] and identify the term 
which contributes to the S-matrix clement under consideration. 

(c) Identify the ten Feynman graphs contributing to the S'-matrix 
element. 
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3. Calculate the symmetry factors for the following graphs: 

(a) 
(b) 
(c) 
(d) 
(e) 
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